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A NEW SCHEME OF INTEGRABILITY FOR 
(BI)HAMILTONIAN PDE 


ALBERTO DE SOLE, VICTOR G. KAC, AND DANIELE VALERI 


Abstract. We develop a new method for constructing integrable Hamilton¬ 
ian hierarchies of Lax type equations, which combines the fractional powers 
technique of Gelfand and Dickey, and the classical Hamiltonian reduction tech¬ 
nique of Drinfeld and Sokolov. The method is based on the notion of an 
Adler type matrix pseudodifferential operator and the notion of a generalized 
quasideterminant. We also introduce the notion of a dispersionless Adler type 
series, which is applied to the study of dispersionless Hamiltonian equations. 
Non-commutative Hamiltonian equations are discussed in this framework as 
well. 
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1. Introduction 

In the present paper we propose a general method of constructing integrable 
(bi)Hamiltonian hierarchies of PDE’s, which combines two most famous approaches. 
The first one is the Gelfand-Dickey technique, based on the Lax pair method, which 
consists in taking fractional powers of pseudodifferential operators [GD76, Dic03], 
and the second one is the classical Hamiltonian reduction technique, combined with 
the Zakharov-Shabat method, as developed by Drinfeld and Sokolov [DS85]. 

The central notion of our method is that of a matrix pseudodifferential operator 
of Adler type , as defined in our paper [DSKV15a], It has been derived there starting 
from Adler’s formula for the second Poisson structure for the M -th KdV hierarchy 
[Adl79] . A pseudodifferential operator A{9) over a differential algebra (V, d) is 
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called of Adler type, with respect to a A-bracket {• a •} on V, if 

{A(;z)aA(u))} = A(w + A + 8)l z {z— w— A — <9) -1 A*(A — z ) 
— A(z)L z (z—w—X—d)~ 1 A(w). 


( 1 . 1 ) 


Recall that a A-bracket is a bilinear (over the base field F) map {■*■}: VxV-> V[A] 
satisfying the sesquilinearity and Leibniz rules axioms (see Section 2.2), A(z) is the 
symbol of A(d), A*(d) is the formal adjoint of A{d), and i z denotes the expansion 
in the geometric series for large z. 

The definition of a rectangular matrix pseudodifferential matrix operator of Adler 
type is similar, see [DSKV15a] or Definition 3.1 of the present paper. 

The first basic property of an Adler type matrix pseudodifferential operator A(d), 
proved in [DSKVISa], is that the A-bracket on V restricted to the subalgebra Vi 
generated by the coefficients of the entries of A(d) satisfies the skewcommutativity 
and Jacobi identity axioms of a Poisson vertex algebra (PVA) (in Section 2.2 we 
recall its definition). Thus an Adler type operator automatically provides Vi with 
a Poisson structure. 

For example, we note in [DSKV15a] that, for each positive integer M, the 
“generic” pseudodifferential operator (resp. “generic” differential operator) 



is of Adler type. Consequently we automatically get on the algebra of differential 
polynomials in {uj \ — oo < j < M} (resp. {uj | 0 < j < M}) the so called 2-nd 
Poisson structure for the KP hierarchy (resp. for the M- th KdV hierarchy). This 
Poisson structure for the M- th KdV hierarchy has been conjectured by Adler in 
[Adl79], and subsequently proved by a lengthy calculation in [GD78]. The second 
Poisson structures for the KP hierarchy were discovered by Radul in [Rad87]. Note 
that these structures are different for different M, though the corresponding hier¬ 
archies of equations are isomorphic. Note also that the 1-st Poisson structure for 
the KP hierarchy was previously discovered in [Wat83], and it is well known to be 
easily derived from the 2-nd structure. In fact, this is the case in general: if A(d) + e 
is of Adler type for any constant e, then we call A(d) of bi-Adler type, and we have 
on Vi a pencil of Poisson structures, see Section 6.1. 

The second basic property of an Adler type M x M-matrix pseudodifferential 
operator A(d ) is that it provides a hierarchy of compatible Lax type equations via 
the method of fractional powers: 



(1.3) 


where B(d) is a fc-th root of A(<9), and k, n are positive integers, see [DSKV15a] and 
Section 5 of the present paper. Here, as usual, the subscript + stands for the positive 
part of a pseudodifferential operator. Of course, in particular, for A(d) = Li(<9) and 
k = 1 (resp. A(d ) = L(m)(S) and k = M), equation (1.3) is Sato’s KP hierarchy 
[Sat81] (resp. Gelfand-Dickey’s M- th KdV hierarchies [GD76]). See also [Dic03]. 

Furthermore, equations (1.3) are Hamiltonian with respect to the Poisson struc¬ 
ture described above, and they are bi-Hamiltonian if A(d) is of bi-Adler type, see 
Section 6. 

The third basic property of Adler type pseudodifferential operators A(d) is that 
it provides an infinite set of conserved densities for the hierarchy (1.3): 



(1.4) 
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Moreover, provided that A(<9) is bi-Adler, these conserved densities satisfy the 
(generalized) Lenard-Magri scheme, introduced in [Mag78]. See [DSKV15a] and 
Sections 5 and 6 of the present paper for details. 

However, our approach allows one not only to recover the basic results of the 
KP theory and M-th KdV theory, including the matrix case, as has been demon¬ 
strated in [DSKV15a], but to go far beyond that in constructing new integrable 
(bi)Hamiltonian hierarchies of Lax type equations. Of course, for this we need to 
construct new Adler type operators. 

Our basic example is the family of N x TV-matrix differential operators 

N 

As (9) = Ijv<9 + 'y ] G Matvrxiv V[9], (1-5) 

i,j =1 

where V is the algebra of differential polynomials in the generators qij , and S t is 
the transpose of a constant matrix S G Mat^xArF. It is not difficult to show, see 
Example 3.4, that the operator As(<9) is of Adler type with respect to the following 
A-bracket on V, making it the affine PVA over the Lie algebra 

{oa&}s = [a, b\ + ti(ab)\ + tr(S 4 [a, &]) , a, b G gl^ . (1.6) 

Unfortunately, the operator As(d) does not have “good” A^-th roots, see Appen¬ 
dix A. The way out is provided by yet another remarkable property of a square 
matrix non-degenerate Adler type operator A(9): its inverse is of Adler type as 
well, with respect to the negative of the A-bracket for A(d) on V. It follows that 
the generalized quasideterminants of A(9) are of Adler type with respect to the 
same A-bracket as for A(9), see Section 4. 

Recall that a quasideterminant of an invertible matrix A over an associative 
(not necessarily commutative) unital ring is the inverse (if it exists) of an entry of 
A -1 . This notion was the key to the systematic development of linear algebra over 
non-commutative associative rings in a series of papers by Gelfand and Retakh, 
and their collaborators, see [GGRW05] for a review. The basic idea of this work 
is that in linear algebra one does not need determinants, all one needs are inverse 
matrices! (We find it surprising that, nevertheless, in the above mentioned review 
the authors of a modern textbook in linear algebra are criticized for the claim that 
“determinants ... are of much less importance that they once were”.) 

In the present paper we use the following generalization of quasideterminants: if 
I G MatjvxMlF and J G MatMxJV®’ are rectangular matrices of rank M < N, the 
(/, J) quasideterminant of an invertible matrix A is |A|/j = (JA -1 /) -1 , assuming 
that the M x M matrix J A -1 1 is invertible. 

To demonstrate our new method of integrability, consider, for example, the ma¬ 
trix A e s defined in (1.5), which is of Adler type with respect to the pencil of PVA 
A-brackets {• a -} e = {• a "lo+ei" a -}i on V, defined by (1.6). Let I G Matjv X A/F and 
J G MatjV/x/v F be rectangular matrices of rank M < N such that S = IJ (/ and J 
are uniquely defined up to a change of basis in F M ). Then, by the results of Section 
4, the generalized quasideterminant |A|/j is a matrix pseudodifferential operator of 
bi-Adler type with respect to the same pencil of PVA A-brackets {• a -} e . Asa conse¬ 
quence, we automatically get, by the above mentioned method of fractional powers 
(cf. (1.4)), a sequence of Hamiltonian functionals fh n = — Resa tr(| A(d)\u) n , 
n > 1, fho = 0, which are in involution with respect to both PVA A-brackets 
{• a *}o and {• a "lii they satisfy the Lenard-Magri recursion relation in the subal¬ 
gebra Vi C V generated by the coefficients of |A|jj, and therefore produce an inte¬ 
grable hierarchy of bi-Hamiltonian equations £jr- = {IbiA M }o| A _ 0 = {h n +l\U}l | a_ 0 ’ 
n > 0, u G Vi. 
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Unfortunately it is a well established tradition to call the Poisson structures 
corresponding to the A-brackets {• a - }o and {• a -}i the second and the first Poisson 
structures respectively, which disagrees with our notation. 

In Section 7 we describe this example in detail in the case of gl 2 , deriving the 
first non-trivial equation of the corresponding hierarchy, for all possible choices of 
the matrix S in (1.5). 

Next, we introduce in Section 8 the notion of a dispersionless Adler type series, 
which allows us to construct integrable Hamiltonian dispersionless hierarchy in a 
unified fashion. For example, the symbols Lm{z) and L^m){z) of the “generic” 
pseudodifferential operators Lm{ d) and L(m)(< 9) given by (1.2) are dispersionless 
Adler type series. The series L\{z) gives rise to the dispersionless KP hierarchy, 
denoted by dKP, introduced in full generality in [KG89], of which the first non¬ 
trivial equation is Benney’s equation for moments, as described by Lebedev and 
Manin in [LM79]. As a result we endow the dKP hierarchy with an infinite series 
of bi-Poisson structures, associated to Lm(z) for every M > 1. The series L(m)(z) 
provides similar results for the dispersionless M-th KdV hierarchy. 

Finally, in Section 9 we develop further the notion of an Adler type pseudodif¬ 
ferential operator with respect to a double PVA, introduced in [DSKV15b], 

As an application of our method, in our next paper [DSKV15new] we will con¬ 
struct an integrable hierarchy of bi-Hamiltonian equations associated to any affine 
W-algebra W(gljv,/)) where / is any nilpotent element of the Lie algebra gl^. 
This was previously known only for nilpotent elements / of very specific type, see 
[DSKV13] and references there. 

Our base held F is a held of characteristic 0. 

The hrst two authors would like to acknowledge the hospitality of IHES, France, 
where this work was completed in the summer of 2015. The third author would 
like to acknowledge the help of the Italian government in the evacuation from the 
earthquake in Nepal in April 2015, as well as the friendship and hospitality of the 
Nepalese people, and the hospitality of the University of Rome La Sapienza during 
his subsequent visit in Rome. 

2. Differential algebras, A-brackets and Poisson vertex algebras 

2.1. Differential algebras and matrix pseudodifferential operators. By a 

differential algebra V we mean a unital commutative associative algebra with a 
derivation d. As usual, we denote by f : V —> V/dV the canonical quotient map. 

Let V((<9 -1 )) be the algebra of scalar pseudodifferential operators over V. Given 
A(d) £ V((<9 -1 )), we denote by A(z) £ V((z -1 )) its symbol, obtained by replacing 
d n (on the right) by z n . The product o on V((d -1 )) can be defined in terms of 
symbols by 

(AoB){z)=A(z + d)B(z). (2.1) 

Here and further, we always expand an expression as ( z+d) n in non-negative powers 
of d. The formal adjoint of the pseudodifferential operator A(d) = '}2^ = _ 00 a n d n 
is dehned as 

N 

A*(d)= J2 

n=— oo 

Introduce an important notation which will be used throughout the paper. Given 
a formal Laurent series A(z) = J2n =-oo a nZ n £ V((z -1 )) and elements b, c £ V, we 
let: 

N 

A ( z + X )(\ x=a b ) c = a n ((z + d) n b)c. 


n=—oo 


( 2 . 2 ) 
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For example, with this notation, the formal adjoint of a pseudodifferential operator 
A(3) £ V((d~ 1 )) has symbol 

A*{z) = {\ x=d A {-z-x)). (2.3) 

Let MatMx a V((9 -1 )) be the space of M x N matrix pseudodifferential op¬ 
erators over V. Let A(d) = (Ay {&)') ,=i . M £ Mat mxJvV(( 3 -1 )) be a matrix 

j = l,...,N 

pseudodifferential operator. Its formal adjoint A*(d) £ MatjvxM V((3 -1 )) has en¬ 
tries ( A*)ij(d) = ( Aji)*(d ). Let I C {1,..., M} and J C {1,..., N} be arbitrary 
subsets. Given an M x N matrix A = (Ay)i=i,...,M, we denote by Ajj the corre¬ 
sponding |/| x |J| submatrix: 


■ 4 « = (Ai), 6W . (2-4) 

Recall that the residue Res z A(z) of a formal Laurent series A(z) £ V((z -1 )) is 
defined as the coefficient of z _1 . We introduce the following notation: 

i z {z- w)~ 1 = ^ z~ n ~ 1 w n , (2.5) 

n£Z_|_ 


i.e. l z denotes the geometric series expansion in the domain of large z (and similarly, 
l w denotes the expansion for large w). According to the previous convention (as in 
(2.1)), a Laurent series in (z + <9) is always expanded using i z . Similarly, if we have 
a Laurent series in (z + A) we omit the symbol l z , and l z or l w will only be used in 
Laurent series involving both variables z and w. Recall that, for a formal Laurent 
series A(z) = J2n=-oo a n zn G V((z -1 )), we have 


N 

Res z A(z)i z (z — w)- 1 = A(w)+ := ^ a n w n . (2.6) 

71 = 0 

We state here some simple facts about matrix pseudodifferential operators which 
will be used later. 

Lemma 2.1. (a) Given two pseudodifferential operators A(d) 1 B(d) £ V((<9 -1 )), 
we have 

Res z A(z)B*( A - z) = Res z A(z + A + d)B(z). (2.7) 

(b) Given two matrix pseudodifferential operators A(d),B(d) £ Mat jv x ^(r 1 )), 
we have 

f Res z tr(A(z + d)B(z)) = J Res z tr (B(z + d)A(z )). (2.8) 

Proof. First, for arbitrary m, n £ Z, we have Res z z m (z — A) n = Res z (z + X) m z n £ 
F[A]((z -1 )) (it reduces to the identity on binomial coefficients = ( — l) fc (fe) 

for k £ Z + ). Hence, by linearity, 

Res z A(z)B(z — A) = Res z A{z + X)B(z). (2.9) 

for every A(z),B(z) £ V((z~ 1 )). On the other hand, in notation (2.2) equation 
(2.7) can be rewritten as 


Res z A(z)(\ x=g B(z - A - x)) = Res z A(z + A + x)(\ x=g B(z)) , 


( 2 . 10 ) 
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which holds by (2.9) (with A + x in place of A). Next, we prove part (b). We have 



= J Res* tr (B(z + d)A(z)). 


In the second equality we used equation (2.7) with A = 0, in the third equality we 
used (2.3), and in the fourth equality we performed integration by parts. □ 

2.2. A-brackets and Poisson vertex algebras. Recall from [BDSK09] that a A- 
bracket on the differential algebra V is a bilinear (over F) map {• a •} : V x V —t V[A], 
satisfying the following axioms (a, b, c £ V): 

(i) sesquilinearity: {da\b} = — A{aA&}, {a\db} = (A + <9){aA&}; 

(ii) Leibniz rules (see notation (2.2)): 

{a\bc} = {a\b}c + {a\c}b , {ab x c} = {aA+xc}(| a;=a &) + {b\ +x c}(\ x=g a,) . 

We say that V is a Poisson vertex algebra (PVA) if the A-bracket {• a •} satisfies 
(a, 6, c G V) 

(hi) skewsymmetry: {b\a} = ~(\ x=g {a-\- x b}); 

(iv) Jacobi identity: {a A {&Ac}} - {^{a ac}} = {{a\b}\ +tJj c}. 

Recall [BDSK09] that if {iqjig/ is a set of generators for the differential algebra V, 
then a A-bracket on V is uniquely determined by {ui\Uj} for i, j £ I. Moreover, if 
the skewsymmetry and Jacobi identity axioms hold on generators, then V is a PVA. 
If a differential algebra V is endowed with a A-bracket satisfying skewsymmetry 
and Jacobi identity, we shall often say that V carries a PVA structure, or simply a 
Poisson structure. 

Example 2.2. Let V be the algebra of differential polynomials over gf v . To keep 
the notation separate, we shall denote by Ejj £ MatArx-/v]F the standard matrix 
with 1 in position (i,j) and 0 everywhere else, and by qy,- £ gf^y the same matrix, 
when viewed as an element of the differential algebra V. Then V = F[<j^ | i, j — 
1,..., N, n £ Z+]. Let S £ MatjvxA F be a fixed matrix and let S t be its transpose. 
The corresponding affine Poisson vertex algebra Vs(g[jv) is the differential algebra 
V endowed with the A-bracket {• a -}s defined by 

{«a^}s = [a, 6] + tr(afc)A + tr(5*[a, b \), a, b £ gl N , (2.11) 

and (uniquely) extended to a PVA A-bracket on V by the sesquilinearity axioms 
and the Leibniz rules. 

2.3. Properties of A-brackets. Later on we shall use the following results from 
[DSK13], which hold on an arbitrary differential algebra V with a A-bracket {• a •}. 
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Lemma 2.3. For A(d),B(d) £ Matjvxiv V((c? 1 )), a £ V, and i, j = 1, ..., TV, we 
have the following identities in V[A](( 2 -1 )) 


N 

{»a(Ao B)ij(z)} = Y^{a\ A ik{z + x)}(\ x=a B kj (z)) 
k=l 
N 

+ A ik(z + A + d){a\Bkj(z)} , 
fe=l 
N 

{(A o B)ij(z) x a} = ^2{B kj (z) x+x a}(\ x=d (A*) ki (X - z)) 
fe=l 
N 

+ ^ ^ {Aj k {z -f- x)x+ x a} (| T _ ftB k j (z)) . 
fc=l 


( 2 . 12 ) 


Proof. It follows immediately from the Leibniz rules and the sesquilinearity axioms 
(cf. the proof of [DSK13, Cor.3.11]). □ 


Corollary 2.4. For A(d) £ MatjvxAf V((9 1 )), a £ V, i,j = 1,... ,7V, andn £ Z+, 
t/ie following identities hold in V[A]((z^ 1 )) 

(a A (A")y(s)} 

n—1 iV 

= E E (^”" 1 "%(* + A + 5){a A A /lfc (* + *)}(| a=8 (^) fca -(^)) ) 

0 h,k=l 

{{A n Uz)xa} 

n —1 iV 

= E E {A fc (^ + x) A+a:+ ,a}(| x=£ ,(A"- 1 -^) fcj -(z))(| 2/=a (^) w (A -*))• 

£=0 h,k—l 

(2.13) 

Proof. It follows immediately from Lemma 2.3. □ 


Lemma 2.5. Let A(<9) £ Matjv X iv V((d -1 )) fee an invertible matrix pseudodiffer¬ 
ential operator, and let A~ l (9) £ Mat at x n V((9 -1 )) be its inverse. Then, for a £ V 
and i,j = 1,... ,7V, the following identities hold in V[A]((,z -1 )) 


{ax(A l )i 3 {z)} 


= ~ E ( A ^)M(z+A+0){a A A, t (,z + a:)}( ^{A , 

S,£=l 
N 

J2{Ast{z+x)x+x+ v a}( ( A~ % ) tj {z )) ( =a (A*- x ) si (A- 2 )) . 


S,t= 1 


Proof. This is the same as [DSK13, Lem.3.9] written in notation (2.2). 


(2.14) 

□ 


Corollary 2.6. Let A[d) £ MatjvxJV V((<9 1 )) be an invertible matrix pseudodif¬ 
ferential operator. Then, for a £ V and i, j = 1,..., TV, and n < —1, the following 
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identities hold inV[X)((z 1 )) 

(a A (A")«(z)} 

-1 N 

= -EE (^- 1 -V(^ + A + 5){a A ^ fc ( 2 + x)}(|_ a (A%(z)), 

£=n h,k =1 

{(A")^) a o} 

-1 TV 

= -EE {^ fc ^ + c C ) A+ , + ,a}(| x=s (A"- 1 ^) fc ,(z))(| y=a (^), l (A-^)). 

^=n h,k= 1 

(2.15) 

Proof. It follows from Corollary 2.4 and Lemma 2.5. □ 


3. Operators of Adler type and their properties 

3.1. Adler type matrix pseudodifferential operators. 

Definition 3.1 (see [DSKV15a]). An M x N matrix pseudodifferential operator 
A{&) over a differential algebra V is of Adler type with respect to a A-bracket {• \ •} 
on V, if, for every (i,j), (h, k) £ {1,..., M} x {1,..., N}, we have 

{Aij(z)\A hk (w)} = A hj (w + A + d)i z (z-w-X-d)~ 1 (A ik )*(X - z) 

i l o • JL ) 

- A hj (z)i z (z-w-X-d)~ A ik (w). 

(In (3.1) (Aik)*(d) denotes the formal adjoint of the scalar pseudodifferential oper¬ 
ator A ik (d), and (A ik )*(z) is its symbol.) 

Note that in [DSKV15a] the Adler type operators are defined with the opposite 
sign in the RHS of (3.1). Obviously, this difference in sign does not affect any of the 
results in that paper, but the present choice of sign seems to be more convenient. 

Remark 3.2. Note that the RHS of (3.1) is in fact regular at z = w + A + d, namely 
we can replace the expansion l z by the expansion i w . To see this, just observe that 

l z (z — w — X — <9) -1 — l w (z — w — A — <9) -1 = S(z — w — A — d ), 

the formal 5-function, and use the following properties of the formal delta function: 


A(z)S(z — w — X — d) = A{w + A + d ), S(z — w — A — 9)A*(A — z) = A{w) 
for any scalar pseudodifferential operator A{d). 

Example 3.3. Any M x N matrix with entries in F is a pseudodifferential operator 
of Adler type for any A-bracket on V. 


Example 3.4. Consider the affine PVA Vs(g( Ar ) from Example 2.2. Consider the 
following N x N matrix differential operator over Vs(flljv): 


N 

As{d) = livd + + S 4 

*> 2=1 


/ d + qn 
qu 

921 

9 + 922 ■ 

■ 9ivi ") 

9N2 

\ <71 N 

q-2N 

d + qNN / 


(3.2) 
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We claim that the matrix differential operator As {&) is of Adler type with respect 
to the affine A-bracket (2.11). Indeed, the RHS of (3.1) is, for the matrix (3.2), 

A hj (w + A + d)L z (z-w-X-d)~ 1 (A ik )*(X - z) - A hj (z)t z (z-w-X- <9) -1 A ik {w) 
= ( 5hj(w + A + d) + qjh + Sjh)i z {z — w — X — d) 1 (Si k {z — A) + q k i + s k i) 

- (S hj z + q jh + Sj h )i z (z - w - X - d^iSikW + q ki + s ki ) 

= Sikiqjh, T Sjh) Shj{q k i T Ski) T A 8 k j8i k 
= {Qji\Qkh}S = {Aij(z)\A hk (w)}s ■ 

Example 3.5. Let M be a positive integer and let V be the algebra of differential 
polynomials in the infinitely many variables Uj- a b , where 1 < a,b < N and —oo < 
j < M (resp. in the finitely many variables Uj ;a b, where 1 < a, b < N and 0 < 
j < M). Consider the following “generic” N x N matrix pseudodifferential (resp. 
differential) operators: 

M M 

Lmn( 9)= ^2 Ujd 3 (resp. L (MN) (d ) = ^Ujd 3 ) with U M = In , 

j=—o o j=0 

(3.3) 

where Uj = \Uj- a b) b _ 1 G Mat^xN V for every j < M. We claim that, for every 
S € MatjvxJvlF, the one parameter family of operators A(<9) = Lmn(9 ) +eS (resp. 
A{9) = L(mn) + oS), e £ F, is of Adler type with respect to a pencil of A-brackets 
{• \ } e = {• a *}o + e{- a -}i on V. Indeed, in this case the LHS of (3.1) becomes a 
Laurent series (resp. polynomial) in z and w of highest degree IV — 1 in z and in 
w, whose coefficients are the unknown A-brackets {ui; a bx u j;cd}e, which are subject 
only to the conditions coming from the skewsymmetry axiom of Section 2.2, and 
to the assumption of being linear in e. On the other hand, the RHS of (3.1) has 
obviously degree at most N — 1 in z, and therefore it has also degree at most N — 1 
in w due to Remark 3.2, it is a quadratic expression in e, and the coefficient of 
e 2 is immediately seen to be zero. In conclusion, the pencil of A-brackets {-A‘}e 
is defined on the generators of V by (3.1), and it extends to the algebra V by the 
sesquilinearity axioms and the Leibniz rules, [BDSK09]. We point out that, as we 
shall show in [DSKV15new], L(m N)(d) is of Adler type with respect to the PVA 
structure of the W-algebra associated to the Lie algebra gl NM and rectangular 
nilpotent / £ gl, VM with N Jordan blocks of size M. 

The relation between Adler type operators and Poisson vertex algebras is given 
by the following result, which is an immediate consequence of [DSKV15a, Lem.4.1] 
(see also [DSKV15a, Rem.2.6]): 

Theorem 3.6. Let A(d) £ MatMxN V((<9 -1 )) be an M x N-matrix pseudodiffer¬ 
ential operator of Adler type with respect to a X-bracket {• a •} on V, and assume 
that the coefficients of the entries of the matrix A(d) generate V as a differential 
algebra. Then V is a Poisson vertex algebra (with X-bracket {• A’}j- 

3.2. Properties of Adler type pseudodifferential operators. 

Proposition 3.7. Let V be a differential algebra with a X-bracket {• a •}• Let A(d) 
be an M x N matrix pseudodifferential operator of Adler type with respect to the 
X-bracket ofV. Then the following properties hold: 

(a) For every subsets I C {1,..., M}, J C {1,..., N }, the corresponding |/| x | J| 
submatrix Ajj is of Adler type (with respect to the same X-bracket {• a *}^- I n 
particular, every entry Aiffd) of the matrix A(d) is a scalar pseudodifferential 
operator of Adler type. 
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(b) The formal adjoint N x M matrix pseudodifferential operator A* (<9) is of Adler 
type (with respect to the same X-bracket {• \ ■}). 

(c) Assuming that M = N and that the square matrix pseudodifferential operator 
A(d) is invertible in the algebra Matjvxiv V((9^ 1 )), the inverse matrix A(d)~ 1 
is of Adler type with respect to the opposite X-bracket —{- ^-}. 

(d) If A(d) G MatMxiv V((<9 -1 )) and B(d) G Matjv x p V((9 -1 )) are both matrices 

of Adler type with respect to the X-bracket {-A - }; and if {Aij(z)\Bh k (w)} = 0 
for all i,j,h,k, then the product A(d) o B(d) G MatMxP is a matrix 

of Adler type. In particular, if S G MatpxM®' andT G MatjvxQlF are constant 
matrices, then the matrix SA(d)T G Matp X Q V((9 -1 )) is of Adler type. 


Proof. Part (a) is an obvious consequence of the definition (3.1) of Adler type matrix 
pseudodifferential operators. Next, let us prove part (b). By the sesquilinearity 
assumptions, we have, using notation (2.2) 

{(A*) i: ,(z)A(A’ ,i )f lfc (w))} = (" {Aji(X - z)\A kh (-w - X - a;)}) . (3.4) 

We can now use condition (3.1) on the matrix A(d) to rewrite the RHS of (3.4) as 

( - A ki {-w - x + d)t z {z - w - X - x + <9) _1 (A jh )*(z) 

+ A ki (X - z)l z (z - w - X - x + d)~ 1 Aj h (—w - A - x)j 

= -(A*) hj (z)i z (z -w-X- d)~ 1 (A*) ik (w) 

+ (. A*) hj (w + X + d)i z (z -w-X- <9) _1 ((A*) ifc )*(A - z ), 

proving that the matrix pseudodifferential operator A* (d) is of Adler type. 

Next, we prove part (c). By Lemma 2.5, we have 

{(A-^MA- 1 )^)} 

N 

— ^ ( (A )/is(lU “t“ X d^jfApqi^Z A X± ) A-f-rci -\-Xs A s t (w A £2)} 

p,q,s,t= 1 

X ( (A-%(*))( (A~ 1 ) tk (w) \ ( (A* _1 ) P i(A - z) \ . 

We then apply condition (3.1) on the matrix A(d) to rewrite the RHS of (3.5) as 

N 


(3.5) 


Y {A 1 )hs(w + X + d) 

p,q,s,t—l 


x ^A sq (w+X+xi+X 2 +x 3 +d)L z (z-w-X-x 2 -x 3 -d) 1 (A pt )*(X+x 3 -z) ^ 
- A sq (z+xi)i z (z-w-X-X 2 -x 3 -d)~ 1 A pt (w + x 2 )j 
X ( (A-y qj (z)) ( (A-Vm) ( (A*-\ z (X - zj) . 

\ X\—0 / \ X-2—O / \ X^—O / 


x\=d 

By definition of inverse matrix we have 
N 

^ 1 )hs{' w + A + 9) o A sq (w + A + d) = Sqh , 


S= 1 

N 


\-\-d—z)(A* 1 ) P i(A — z) = Sa , 
p =1 

N N 

Y A sq( z +d)(A~ 1 ) qj (z) = S sj , Y A pt( w + d)(A~ 1 ) tk {w) = 8pk . 


g=l 
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Hence, the RHS of (3.6) can be rewritten as 

(A~ 1 ) hj (z)i z (z - w - A - <9) _1 (A _1 ) ife (w) 

- (A~ 1 ) hj (w + A + d)i x (z — w — X — d)~ 1 (A*~ 1 ) ki (X - z ), 

proving that A _1 (<9) is of Adler type with respect to the opposite A-bracket —{- a •}- 
Finally, the proof of (d) is analogous to the proof of [DSKV15a, Prop.2.14]. We 
reproduce here the argument for completeness. By Lemma 2.3 and the assumption 
that the A-bracket between entries of A and B is zero, we have 
{(A o B)ij(z)\(A o B) hk (w)} 

N 

= i A is( z + x) x+x A ht (w + y)}(\ x=a Bsj{z))(\ y=a B tk {w)) 

s,t= i (3.7) 

N 

+ Aht{w + A + d){B S j(z)\+ x B t k(w)}(\ x _ a (Ai S )*(\ — z)) . 

S,t= 1 

By the Adler condition (3.1), the first term in the RHS of (3.7) is equal to 
N 

y. A hs (w + A + d)B sj (z)t z (z -w - A - S) _1 (A it )*(A - z)B tk (w) 

«,t=i ( 3 - 8 ) 

- (Ao B) hj (z)t z (z - w - A - 9) _1 (Ao B) ik {w ), 
while the second term in the RHS of (3.7) is 

(A o B) hj (w + A + d)i z (z — w — X — 9) _1 ((A o B) ik )*{X - z) 

N (3.9) 

- y A ht (w + A + d)B t j(z)i~(z -w-X - d)~ 1 (A is )*(X - z)B sk (w). 

S,t=l 

Combining (3.8) and (3.9), we conclude that Ao B satisfies the Adler condition 
(3.1), as claimed. □ 

4. Quasideterminants of operators of Adler type 

4.1. A brief theory of quasideterminants. In this section we review the notion 
of quasideterminants of a matrix over an arbitrary unital associative ring, and 
we prove some of their properties. For an extended review on the subject, see 
[GGRW05|. 

Let R be a unital associative ring. Let A £ Matjv x ./v R -be an N x N matrix with 
entries in R. Let I, J C {1,..., N} be subsets of the same cardinality |/| = | J\ = M. 

Definition 4.1. The (/, J)-quasideterminant of A is defined as (recall the notation 
(2.4)) 

\A\ij'.= ((A : )j/) 1 £ Mat mxmR, (4.1) 

assuming that the RHS makes sense, i.e. that the N x N matrix A is invertible, 
and that the M x M matrix (A -1 )jj is also invertible. (The Gelfand-Retakh 
quasideterminant was originally defined in the special case M = 1.) 

When the quasideterminant is defined, it is possible to write a simple formula 
for it. This is given by the following 

Proposition 4.2. Assume that A £ Mat^xiv R is invertible and that the submatrix 
A/cjc £ Mat(jv_M)x(A'-M) Is invertible, where I c and J c denote the complement 
of the sets I and J in {1,, N }. Then the (J, J)-quasideterminant of A exists, 
and it is given by the following formula: 

\A\u = A/j — A IJ c(A I cjc)~ 1 A I cj . 


(4.2) 
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Proof. Without loss of generality, we can assume that I = J = {1 C 

{1,,N}. In this case, the matrix A can be written in blocks as follows: 


A = 


Ajj 

A/cj 


Ajj? 

Ajcjc 


Accordingly, the inverse matrix has the block form 


a-i _ ( {A~ 1 ) JI (A~ 1 ) JI c 

- V ( A(A- 1 )^ 

Under the assumptions that A and Ajcjc are invertible, it is easy to find a formula 
for the four blocks of A -1 : 

(A~ 1 )ji = (Ajj — Aijc(Ajcjc)~ 1 Aicj )~ 1 , 

( A~ 1 ) JI c = - (A^ 1 )jjAjjc(Ajcjc )~ 1 , 

(A~ 1 )jc I = —(Aicjc)~ 1 Aicj(A^ 1 )jj , 

(A- 1 ) jo jo = (Ajojc )- 1 + ( Aicjc)- 1 Ajcj(A~ 1 )jiAjjc(Ajcjc ) _1 . 

The claim follows. □ 



Quasideterminant satisfy the so-called hereditary property , which can be stated 
as follows: 

Proposition 4.3. Let I\ C I C {1,..., TV} and J\ C J C {1,..., TV} be subsets 
such that |/i| = |Ji| and |/| = |J|. Assuming that the (I , J )-quasideterminant and 
the (I i, Ji )-quasideterminant of A exist, we have 

\A\i 1 j 1 = \\A\ij\i iJi . ( 4 . 3 ) 

Proof. By the definition (4.1) of quasideterminant we have 

\\ A \ ij \ i iJi = ((((( A - 1 )^)” 1 )- 1 )^)- 1 = ((A- 1 )**)" 1 = |A| 7lJl . 

□ 


4.2. Generalized quasideterminants. We can generalize the notion of quaside¬ 
terminant as follows. Let, as before, A £ Mat nxn R, and let I £ Mat nxm 
J £ MatMxA -R, for some M < N. 

Definition 4.4. The (/, J)-quasideterminant of A is 

\A\u = (JA" 1 /)- 1 £ Mat mxm R , (4.4) 

assuming that the RHS makes sense, i.e. that A is invertible in Mat^xArR and 
that JA -1 / is invertible in Matjuxif R- 

Note that Definition 4.4 is a generalization of Definition 4.1. Indeed, given the 
subsets / = {ii, ■ ■ ■ Am} and J = {j i,... , Jm} of {1,. .., N}, then the quasideter¬ 
minant \A\u defined by (4.1) coincides with the quasideterminant \A\c^Rj defined 
by (4.4), where Cj is the TV x M matrix with 1 in position (ie,£), for £ = 1,..., M, 
and zero otherwise, and Rj is the M x TV matrix with 1 in position (l, jj), for 
£= 1,..., M, and zero otherwise. The following result, though very simple, will be 
very useful in our theory. 

Theorem 4.5. Let A £ MativxivR, I G Matj\ r X M R, J £ Mat mxn R, for some 
M < TV. Assume that the (I , J)- quasideterminant |A|/j exists and that the matrix 
A + IJ £ MatTVxivR is invertible. Then, the (I , J)-quasideterminant of A + IJ 
exists, and it is given by 


IA + IJ\ij = \A\u + • 


(4.5) 
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Proof. If we multiply the identity A +1J = A +1J by J (A +1 J) 1 on the left and 
by A~ 1 I\A\ij on the right, we get 

J{A + IJ) 1 /(1m + |A|/j) = 1m , 

i.e. 1 m + \A\iJ i s a right inverse of J(A + IJ)~ 1 I. Similarly, if we multiply the 
same identity A + IJ = A + IJ by \A\uJA~ X on the left and by (A + IJ ) -1 / on 
the right, we get 

(1m + \A\ij)J(A + IJ) 1 J = Ijvf , 

i.e. 1 m + \A\iJ i s a l so a left inverse of J(A + /J) -1 /. In conclusion, 1 m + |A|/j = 
(J(A + IJ)~ 1 I)~ 1 = \A + IJ\ij, proving the claim. □ 

4.3. Adler type operators and quasideterminants. 

Proposition 4.6. Let V be a differential algebra with a X-bracket {•>•}. Let A(d) £ 
MatjvxA V((<9 -1 )) be a matrix pseudodifferential operator of Adler type with respect 
to the X-bracket ofV. Then, for every I £ MatjvxMlF and J £ MatMxAr®’ with 
M < N, the quasideterminant \A(d)\u is a matrix pseudodifferential operator of 
Adler type, provided that it exists. 

Proof. It is an obvious consequence of Proposition 3.7(c) and (d) and of the defini¬ 
tion (4.4) of quasideterminant. □ 

Remark 4.7. For I, J C {l,...,iV} of the same cardinality, one can show, by 

direct computations similar to the ones used to prove Proposition 3.7(c), that 
Aij(d) — Aijc(d) o (Ajcjc(<9)) -1 o Aicj(d) is a matrix of Alder type, even if the 
quasideterminant \A(d)\u does not exist. We omit the proof of this fact. 

5. Integrable hierarchy associated to a matrix pseudodifferential 

OPERATOR OF ADLER TYPE 

Let V be a Poisson vertex algebra with A-bracket {•>•}■ We have the corre¬ 
sponding Lie algebra structure on V /dV with Lie bracket 

Uf>fff} = f{hff} | A =o> (5- 1 ) 

and a representation of the Lie algebra V /dV on V given by the action 

{//>5} = {/a 5}| a=0 - (5-2) 

Recall that the basic problem in the theory of integrability is to construct an infinite 
linearly independent sequence of elements Jh n £ V/dV, n £ Z+, called Hamiltonian 
functionals, which are in involution, i.e. such that 

{fhm,fhn} = 0 for all m, n £ Z+ . 

In this case we obtain a hierarchy of compatible Hamiltonian equations 

= {fK,u}, u £ V. 

In the present section we show that, given a matrix pseudodifferential operator 
of Adler type, it is possible to construct a sequence of Hamiltonian functionals in 
involution (the question of their linearly independence is treated separately). This 
is stated in the following: 

Theorem 5.1. Let V be a differential algebra with a X-bracket {■ \ •}. Let A(<9) £ 
MatjvxA V((<9 -1 )) be a matrix pseudodifferential operator of Adler type with respect 
to the X-bracket {• A •}, and assume that A(d) is invertible in Matjvx nV((5” 1 )). 
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For B(d) £ MatjvxivV((9 1 )) a K-th root of A (i.e. A(d) = B(d) K for K £ 
7\{0}) define the elements h n ,B € V, n £ Z, by 


h n ,B = -j—r Res z tT(B n (z)) forn^O, h 0 = 0 . 
\n\ 


Then: 

(a) All the elements Jh Ut B are Hamiltonian functionals in involution: 


{/ h m ,B, f h Ut c} = 0 for all m , n £ Z, B,C roots of A . 


(5.3) 


(5.4) 


(b) The corresponding compatible hierarchy of Hamiltonian equations satisfies 

“ 77 — = {fh n ,B,A(z)} = [(B n ) + ,A\(z), n £ Z, B root of A (5.5) 
atn,B 

(in the RHS we are taking the symbol of the commutator of matrix pseudodif¬ 
ferential operators), and the Hamiltonian functionals fh n ,c> fi £ Z +; C root of 
A, are integrals of motion of all these equations. 


Remark 5.2. Since, by assumption, A(d) is invertible in Mat^vxAf "^((S -1 )), then, 
for K > 1, a Ji-th root B(d) of A(d) is invertible as well, its inverse being i? -1 = 
A- 1 o B I< ~ 1 . Hence, (5.3) makes sense for all n £ Z. On the other hand, the 
sequence (5.3) for B is related to the same sequence for £? -1 by the following 
simple relation: = h n b -i for all n £ Z+. 


Remark 5.3. One can state the same Theorem 5.1 without the assumption that 
A (and therefore B) is invertible, at the price of assuming that K > 1, and of 
restricting the sequence h n in (5.3) to n £ Z+. Moreover, since the proof of (5.4) 
and (5.5) is based on Lemma 5.5, one needs to restrict equation (5.4) to m > K 
and n > L, and equation (5.5) to n > K. 


Remark 5.4. Let Vi C V be the differential subalgebra generated by the coeffi¬ 
cients of A(d). Then by Theorem 3.6 {• a - } is a PVA A-bracket on Vi, and (5.5) 
uniquely extends to a hierarchy of compatible Hamiltonian equations over the PVA 
Vi. However, if the A-bracket on V is a PVA A-bracket, then we have a hierarchy 
of compatible Hamiltonian equations with Hamiltonian functionals fh Hi B, n £ Z, 
B a A-th root of A, on the whole V. 


In the remainder of the section we will give a proof of Theorem 5.1. It is based 
on the following: 

Lemma 5.5. For a £ V, n £ Z, and Bid) £ Mat^x^ V((i9 -1 )) a K-th root of 
A(d) (K eZ\{0}j, we have 

N 

{^n,B a^I | A=0 = ~ ^ ' R- es z{Aij(z + x) x a} )ji { z )) > 

1 , 7=1 

„ ( 5 - 6 > 
/{<uVb}|;, = 0 = - Y2 I Res *{‘ 1 *Aij(z + x)}| A=0 (|^ a {B n . 

i,j =1 j 

Proof The proof of the first equation in (5.6) is similar to the proof of [DSKV15a, 
Lem.3.2]. We review here the argument. By the second equation in (2.13), we have, 
for n > 0, 

{(B n Uz) x a} 

n-i N / r _n 

r=0 h,k =1 
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while, by the second equation in (2.15), we have, for n < — 1, 

{{B n Uz) x a} 

-1 N 

-EE {^(^ + ^A + x + W(| x= a( s "" 1 “')^W)(L= 9 (^)«(A-^)). 


£=nh,k=l 


(5.8) 

For K > 1, we take n = K in (5.7), replace z by z -\- d and A by d acting on 
B r ]'- K (z), take residue in z and sum over i,j = l,...,N. As a result, we get 


N 

£ Res 2 {Ay(z + X) x a}(\ x=d {B n - K )ji(z)) 
i,j =1 

K- 1 N 

= £ £ 

^=0 i,j,h,k— 1 

X (| x=a (^- 1 ^) fej (z + A))(| y=a (^) w (-^))(| A= 9 J B^(^)) 

K- 1 N 

= £ £ ReSjjB^^^z + A + a: + y)A+x+ya} 

■£=0 Z,J,/i,fe=l 

x (| x=d {B K ~^) kj {z + A + y))(\^ d Bl~ K {z + y)) {\ d (B e ) ih (z)) 


(5.9) 


I x=d '• 
N 


= K Y Res z {B hk (z + x) x a}(\ x=g (B n 1 ) k h(z)) ■ 


h,k=l 

In the second equality we used Lemma 2.1(a). For K < —1, we do the same 
manipulation starting from equation (5.8) in place of (5.7). As a result we get the 
same final equation as (5.9) (which thus holds for every K £ Z\{0}): 

N 

£ Res z {A i:j {z + A) A a}(| x=g (B n ~ K )ji(z)) 


i,j =1 


N 


(5.10) 


= I< Y Res z {B hk (z + x) x a}(\ x=g (B n ^khiz)) ■ 

h,k=l 

On the other hand, for n > 1 we put A = 0 in (5.7), take the residue in z, let i = j 
and sum over i = 1,... ,N, and as a result we get, by the definition (5.3) of 

-K N 


{hn,B x a}\ x=0 =-^Res 2 {(B") ii ( 2 r) A a}| A=0 


—K 

n 


-K 

n 


i =1 

n —1 N 

£ £ Res 2 {B ?lfc ( 2 + x) ;c+y a}(| ;c=a (^- 1 ^) fci (z))(| 2/=9 (B^)^(--)) 

£=0 i,h,k= 1 

n— 1 N 


-£ £ Res z {B hk (z+x+y) x+y a}(\ x=g {B n 1 e ) ki (z + y)) (\ y=d {B% h (z)) 
£—0i : h,k=l 
N 

=-K Y Res z {B hk (z + x) x a}(\ x=g {B n ~ 1 ) kh (z)) . 


h,k= 1 


(5.11) 


Again, in the third equality we used Lemma 2.1(a). We arrive at the same con¬ 
clusion, for n < —1, using equation (5.8) in place of (5.7). Combining equations 
(5.10) and (5.11), we get the first equation in (5.6). 
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Next, we prove the second equation in (5.6). By the first equation in (2.13), we 
have, for n > 1, 

n— 1 N 

{a\{B n )ij(z)} = E E C B n - l ~% h {z + A + d){a x B hk {z + z)}(| x=9 (R%(z)) , 

£—0 h,k=l 

(5.12) 

while by the first equation in (2.15), we have, for n < — 1, 

-1 N 

{a x (B n ) ij (z)} = -J2 E (B n - l -%(z + \ + d){a x B hk (z + x)}(l =a (B%(z)). 

£—n h,k— 1 

_(5.13) 

For K > 1, we let n = K in (5.12), let A = 0, replace zby z+d acting on 71” K ( 2 ), 
take residue in z, apply /, and sum over i,j = l,...,N. As a result we get 

N 


E /Res z {a A (^)^(z + x)}| A=0 (| x=9J B^(z)) 

K -1 AT 

= E E / Re s z (-B Jtf " 1 ^) i?l (z + a){aAB /jfc (z + a;)}| A=0 (| x=a (B n+i -' R: ) fci (z)) 

K-l AT „ 

= E E / R es z (71”+ £ -^) fei (z + 9)(il K - 1 -^(z + 9){a A ^ fc (z)}| A=0 

^=0 i,h,k=l J 
N f 

= I<J2 Res z (il" _1 )fcii(z + 0){a A fl /lfc (z)}| A5=o . 

z, 7 _ 1 • 


A,fc=l * 


(5.14) 

In the second equality we used Lemma 2.1(b). For K < —1, we do the same 
manipulations with equation (5.13) in place of (5.12) and, as a result, we get the 
same final equation as (5.14) (which thus holds for evey K £ Z\{0}): 


N r - 

E J Res z {o A (B if )y(z + a;)}| A=0 (| x=a B}‘ K (z)) 


i,j =1 


N r 

= R E / Res ziB™- 1 )^ + d){a x B hk (z)}\ x=0 

u _ 1 J 


(5.15) 


h.k—1 


On the other hand, for n > 1 we put A = 0 in (5.12), take the residue in z, apply 
J, let i = j and sum over 7 = 1 ,..., TV. As a result we get, by the definition (5.3) 

of h n ^B 1 

—K -2 . f 

f{a X h n ,B }| a= q = 'y ] / Res z {a A (71 ) m( 2 )}| a= q 

n i=i J 

jy. n—1 N „ 

= — E E Res ^ Bn ~ 1 ~ e )ih(z + d){a x B hk {z + x)}\ x=0 (\ x=g {B e ) ki {z)) 

£=0 i,h,k =1 
jr n—1 N ,, 

= yE E / Re s z{B t ) ki { Z + d^-^Uiz + d){a x B hk (z)}\ x=0 

£—0 i,h,k=1 

N r 

= -K ^2 / Res 2 (R” _1 )fc/ t (z + d){a x B hk (z)}\ x ^ 0 . 

u i„ _ 1 J 


h.k— 1 * 


(5.16) 
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Again, in the third equality we used Lemma 2.1(b). For n < — 1 we arrive at the 
same conclusion, using equation (5.13) in place of (5.12). Combining equations 
(5.15) and (5.16), we get the second equation in (5.6). □ 

Proof of Theorem 5.1. Suppose B is a A'-th root of A , AT € Z\{0} and C is an L-th 
root of A, L £ Z\{0}. Applying the second equation in (5.6) first, and then the 
first equation in (5.6), we get 

N r 

^ ^ / Res z \ Ajj(z T x) x A kk (w T y)} 

i,j,h,k=i J ( 5 - 17 ) 

We can now use the Adler condition (3.1) to rewrite the RHS of (5.17) as ciao 

N r 

Res z Res w A h j(w+d)(B m K )ji(z)i z (z-w-d) 1 {A ik )*(-z)(C n L ) k h{w) 

i,j,h,k=l 

N r 

- y Res z Res w (B m ) hi (z)i z (z-w-d) 1 {C n ) ih [w). 

i,h=l ^ 

(5.18) 


For the second term we used the fact that A = B K = C L . By (2.6), we have 
N 

yRes z (B m ~ K )ji(z)L z {z -w- dy 1 o ( A ik )*(-z) = ( B m ) jk (w + d)+ . (5.19) 

i=1 

Hence, the first term of (5.18) is equal to 

N r 

y / Res w A hj {w + d){B m ) jk {w + d)+(C n ~ L ) kh {w) 

j,h,k =1 J 
N 

= y J Res w (C n ) kj (w + d)(B m ) jk (w)+. 


(5.20) 


j,k-l 1 


Here we used Lemma 2.1(b) and the identity C L = A. On the other hand, by (2.6) 
the second term of (5.18) is equal to 

N r 

-y Res w (B m ) hl (w + d)+(C n ) ih (w) (5.21) 

i,h=l ^ 


which, by Lemma 2.1(b), is the same as the RHS of (5.20), with opposite sign. In 
conclusion, (5.18) is zero, as claimed. This proves part (a). 

For part (b), by the first equation in Lemma 5.6 we have 

{fhn jB ,A hk (w)} = {h n ,B x A hk (w )}\ x=0 
N 

= ~y R^ z {Aij(z + x) x A hk (w)}{\ x ^ d (B n - K ) ji (z)) . 

*.j=i 


(5.22) 
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By the Adler condition (3.1) the RHS of (5.22) becomes 


N 


- y; Res z A hj (w + d)(B n K )ji{z)i z {z-w-d) 1 (A ik )*(-z) 

i,3 =1 

N 

+ 22 Res 2 ( A hj{z + d){B n ~ K ) ji (z))L z (z-w-d)~ 1 A ik (w) 
i,j =1 
N 

= -'22 ReSz A hj( w + d)( y B n ) jk {z)i z (z-w)~ l 
i=i 

N 

+ 22 Res z (B n ) hi (z)t z {z-w-d)~ 1 A ik (w) 


(5.23) 


i= 1 


N 

= 22 ( - + <9)(R")ifc(w)+ + (5")w(w + <9) + A ifc (w;)) . 

2=1 

In the second equality we used Lemma 2.1 and the identity A = B K , while in 
the last equality we used (2.6). This proves (5.5) and completes the proof of the 
Theorem. □ 

6. Operators of bi-Adler type and bi-Hamiltonian hierarchies 

6.1. Operators of bi-Adler type and bi-Poisson vertex algebras. Let {• a -}o 

and {• > -}i be two A-brackets on the same differential algebra V. We can consider 
the pencil of A-brackets 

{• a‘} e = {• a - } o + e{‘ a - } i , e€F. (6-1) 

We say that V is a bi-PVA if {• a -} e is a PVA A-bracket on V for every e G F. Clearly, 
for this it suffices that {• a -jo, {• a -}i and {• a -}o + {• a -}i are PVA A-brackets. 

Definition 6.1. Let A(d) G MatMxA V((9 -1 )) be a matrix pseudodifferential 
operator, and let S G MatMxivF. We say that A is of S-Adler type with respect 
to the two A-brackets {• a -}o and {• a -}i if, for every e G F, A(d) + eS is a matrix 
of Adler type with respect to the A-bracket {• a-} e - This is equivalent to saying 
that A(d) is a matrix of Adler type with respect to the A-bracket {• a -}o) he. (3.1) 
holds, and that 

{A ij {z)\A hk (w)} 1 = l z (z-w-\)~ 1 S ik (A hj (w + A) - A hj (z)) 

+ i z (z-w-X-d)~ 1 S hj ((A ik )*(X - z) - A ik (w )) . 

In the case M = N, we also say that A is of bi-Adler type if it is of Iv-Adler type. 
In this case condition (6.2) becomes: 

{Aij(z)\A hk (w)}i = 6 ik i z (z-w-X)~ 1 (A hj (w + A) - A hj (z )) 

+ 6 hj L z (z-w-X-d)~ 1 ((A ik )*(X - z) - A ik {w )) . 

Example 6.2. From Example 3.4 we have that, for every S G MatjvxAfF, the 
matrix differential operator A(d) = 1 jyd + Ei J= i i s °f 5-Adler type with 

respect the the bi-PVA A-brackets (cf. (2.11)) 

{a\b} 0 = [a,b]+tr(ab)X , {aA&}i = tr(S[a, 6]) , a, & G fl l N , (6.4) 

extended to V^gljy) by sesquilinearity and the Leibniz rules. 

Operators of S'-Adler type are related to bi-PVA’s in the same way as operators 
of Adler type are related to PVA’s. In other words, we have the following extension 
of Theorem 3.6, which, in fact, is a consequence of Theorem 3.6: 


( 6 . 2 ) 


(6.3) 
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Theorem 6.3. Let S £ Mat«-xJV F, and let A{d) £ Mat mxn V((d~ 1 )) be an 
M x N-matrix pseudodifferential operator of S-Adler type with respect to the X- 
brackets {• a - } o an d {• a *}i on V. Assume that the coefficients of the entries of 
the matrix A(d) generate V as a differential algebra. Then V is a bi-PVA with the 
X-brackets {• a -}o an d {■ \ -}i. 


6.2. Bi-Hamiltonian hierarchy associated to an operator of bi-Adler type. 

Let V be a bi-Poisson vertex algebra with A-brackets {• a "}o and {• a "}i• A bi- 
Hamiltonian equation is an evolution equation which can be written in Hamilton¬ 
ian form with respect to both PVA A-brackets and two Hamiltonian functionals 
fh 0 ,fh 1 £V/dV: 


du 

dt 


= {Jho,u } 0 = , 


u£V. 


The usual way to prove integrability for a bi-Hamiltonian equation is to solve the 
so called Lenard-Magri recurrence relation (u £ V): 


{fh n ,u}o = {fh n+ 1 ,u} 1 , n£ Z+. 


(6.5) 


In this way, we get the corresponding hierarchy of bi-Hamiltonian equations 
du 

-T— = {fhn, uj o = {Jh n+ I,u}i , n £ Z+, u £ V. 

CLtn 

We next show that, given an operator A(d) £ Mat nxn V((c> _1 )) of bi-Adler type, 
the sequence of local functionals provided by Theorem 5.1 satisfies a generalized 
Lenard-Magri recurrence relation, and therefore gives a hierarchy of bi-Hamiltonian 
equations. 


Theorem 6.4. Let V be a differential algebra with two X-brackets {• a -}o an d 
{• a -}i. Let A(d) £ MatjvxAf V((<9 -1 )) be a matrix pseudodifferential operator of bi- 
Adler type with respect to the X-brackets {• a -}o and {• a -}i, and assume that A(d) 
is invertible in Matjvx a rV((9 1 )). Let B(d) £ MatAxivV((id 1 )) be a K-t.h root 
of A. Then, the elements h n ,B £ V, n £ Z + , given by (5.3) satisfy the following 
generalized Lenard-Magri recurrence relation: 

{fhn t B, A(z)}o = {Jh n+ K,B, A{z)}x = [(£")+, A](z) , n£Z. (6.6) 

Hence, (5.5) is a compatible hierarchy of bi-Hamiltonian equations on the bi-PVA V± 
generated by the coefficients of the entries of A(z). Moreover, all the Hamiltonian 
functionals fh n ,c> n £ Z+, C root of A, are integrals of motion of all the equations 
of this hierarchy. 


Remark 6.5. As in Remark 5.3, one can state the same Theorem 6.4 without the 
assumption that A is invertible. In this case we need to restrict to I\ > 1, and the 
proof of the Lenard-Magri recurrence relation (6.6) only works for n > K. 

Remark 6.6. By Theorem 3.6 the A-brackets {• a - } o an d {- a *}i define a bi-PVA 
structure on Vi. The Lenard-Magri recurrence {Jh n ^B,u }o = {fh n +K,B, u}i holds 
for every n £ Z + and u £ Vi, and (5.5) uniquely extends to a hierarchy of com¬ 
patible bi-Hamiltonian equations over the bi-PVA Vi- However, the Lenard-Magri 
recurrence equation does not necessarily hold for every u £ V even if the two A- 
brackets {• a -}o an d {- a *}i of V are compatible PVA A-brackets. This phenomenon 
does actually occur even in the simplest example discussed in Section 7.1. 

Remark 6.7. Suppose that V is a differential field and A(d) £ MatjvxA r V((<9 -1 )) is 
an operator of bi-Adler type with respect to two A-brackets {• a }o an d {• a }i- Then 
A{9) + eljv is again an operator of bi-Adler type, with respect to the A-brackets 
{• a }o + e{- a }i and {• a }i- Note that A(d) + el a is invertible for all but finitely 
many values of e £ F. Hence, in Theorem 6.4 we do not need to require that A(d) 
is invertible. 
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Proof of Theorem 6.f. By Theorem 5.1(b), we have 

{fh n , B ,A hk (w )}o = [(B n ) + , A} hk (w) , n £ Z+ . 

On the other hand, by the first equation in (5.6) we have 

{ — \hn,B \A-hk (^)} 1 1 q 

N 

= - ReSz{A ij (z+x) x A hk (w)}i(\ x=g (B n ~ K ) ji (z)) . 
i,j =1 


Applying the bi-Adler identity (6.3) we can rewrite the RHS of (6.8) as 

N 

- Resz t z (z-w)~ 1 (A hj (w + d)~ A hj {z + d))(B n ~ K ) jk (z) 

3 =1 
N 

- ^Res z (B n ~ K ) hi (z)t z (z-w-dy 1 ((A ik )*(-z) - A lk (w)) 

i=1 

N N 

= ~y~] A^-(u) + d)(B n ~ K ) jk (w ) + ^](B"“ a )h l (w + <9)+A lfc (w;) 

j=i j=i 

= [(R"- if ) + ,A] hfe H. 


(6.7) 


( 6 . 8 ) 


(6.9) 


In the first equality we used equations (2.6) and (5.19) and the identities AoB n ~ K = 
B n = B n ~ K o A. Comparing (6.7) and (6.9), we get (6.6). The last statement of 
the theorem is an immediate consequence of Theorem 5.1 and (6.6). □ 


6.3. A new scheme of integrability of bi-Hamiltonian PDE. Based on The¬ 
orems 4.5, 5.1 and 6.4, we propose the following scheme of integrability. Let V 
be a differential algebra with compatible PVA A-brackets {• a ’}o and {- a ‘}i- Let 
S £ MatATxArF and let A(d ) £ Matjvx n V((9 -1 )) be an operator of S-Adler type 
with respect to the A-brackets {• a "}o and {- a *}i• Assume (without loss of gener¬ 
ality) that the differential algebra V is generated by the coefficients of A(d). Then, 
we obtain an integrable hierarchy of bi-Hamiltonian equations as follows: 

1. consider the canonical factorization S = IJ, where J : -» Im(S') and I : 

Im(5) ^ F^; 

2. assume that the (/, J)-quasideterminant \A\jj(d) exists; then, by Theorem 4.5 
and Proposition 4.6 | A\jj is an M x M matrix pseudodifferential operator (where 
M = dimlm(S')) of bi-Adler type with respect to the same A-brackets {• A’}o 
and {• a -}i; 

3. consider the family of local functionals {fh n B \n £ Z, B a A'-th root of A} 
given by (5.3); then, by Theorem 6.4, they are all Hamiltonian functionals in 
involution with respect to both PVA A-brackets {• a-} o an d {- a ‘}i? an d they 
satisfy the Lenard-Magri recurrence relation (6.6); 

4. we thus get an integrable hierarchy of bi-Hamiltonian equations 

du 

-j. -= {Jh ni B,u }o = {fh n+ K,B,u }i j (6.10) 

dtn,B 

provided that the Jh Ut B span an infinite dimensional space. 

In Section 7 we will show how to implement this scheme to the affine bi-PVA’s for 
gljY defined in Example 2.2 and the corresponding operators of bi-Adler type in 
Example 3.4. 

In forthcoming publications we plan to implement the same scheme to construct 
integrable hierarchies associated to the classical W-algebras W(g,/) for classical 
Lie algebras g and arbitrary nilpotent elements / £ g. 
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Remark 6 . 8 . Consider the space C\ of Casimir functionals with respect to the A- 
bracket {• a -}i: 

C' 1 = {//ieV/av|{/h,V}i = o}. (6.11) 

Note that, if A{d) is an operator of bi-Adler type, then fhi : A = /Res z trA(,z) 
is a Casimir functional (this follows by Theorem 6.4, or it can be easily checked 
directly). By Theorem 6.4, for this Casimir element the Lenard-Magri recurrence 
scheme (6.5) is infinite, namely it admits an infinite sequence {Jh n } ne i + of solu¬ 
tions, starting from fh o = JhiA- We conjecture that, if A(<9) is a matrix pseudo¬ 
differential operator of bi-Adler type, then the Lenard-Magri scheme (6.5) is infinite 
for each Casimir functional fho £ C\. Note that, by [BDSK09, Prop.2.10(c)] all 
the resulting Hamiltonian functionals are in involution. 

Remark 6.9. Sometimes A(d) is of bi-Adler type with respect to a certain bi-PVA 
V, and the differential subalgebra Vi C V generated by its coefficients is a proper 
bi-PVA subalgebra. In such case it may happen, as the example discussed in Section 
7.1 shows, that fhi.A = /Res z trA(z) is a Casimir functional with respect to the 
A-bracket {• a ’}i on Vi but not on V. 


7. Example: bi-Hamiltonian hierarchies for the affine PVA Vs(gl 2 ) 

Consider the affine PVA V = Vs(gl 2 ) from Example 2.2 and the matrix differ¬ 
ential operator of Adler type A(<9) + S' 4 , where A(d) = d\ 2 + ( ^ 11 ® 21 ), from 

V 912 922 ) 

Example 3.4. We shall study the bi-Hamiltonian hierarchies associated to this op¬ 
erator of S-Adler type, for various choices of S € Mat 2 X 2 F, following the general 
method discussed in Section 6 . 

Note that if S is non-degenerate, then its canonical decomposition is S = Sl 2 
and the corresponding (Sl 2 )-quasideterminant of A(d) is |A(<9)|si 2 = S' 
which is a differential operator of order 1 . Hence, the residue of any power of it is 
zero. Also, by the results of Appendix A, taking square roots of A(d) only seem to 
lead to trivial hierarchies as well (cf. Remark A.5). 

Hence, we only need to consider the case when S has rank 1 . In Section 7.1 
we consider the case when S is diagonalizable of rank 1 , while in Section 7.2 we 
consider the case when S is nilpotent of rank 1 . 


7.1. First case: S = diag(s, 0), s £ F\{0}. We know that the family of Adler type 
operators A(d) + eS, e £ F\{0}, defines two compatible PVA A-brackets {• a -}o and 
{• a -}i on V, given by (cf. (2.11)) 


{(Jij\Qhk\ 0 — 3jhQik $kiQhj T ^jh^ik^ 1 3 L ji h,k £ {T, 2 } , 
and the only non-zero A-brackets among generators for {• a - }i are 


(7.1) 


{ 9 i 2 a 92 i}i — s , { 921 a 9 i 2 jt — — s . (7.2) 

Let us consider the canonical factorization B = IJ where I = ^ ^ ^ and J = 
( s 0 ), and define 

L(d) = |A(<9)|/J = i|A(<9)|ii = ^ (d + q n - q 2 i{d + q 2 2 )~ 1 qi 2 ) ■ (7-3) 

Let Vi C V be the differential subalgebra generated by the coefficients of L{d). 
By Thorem 6.3 it is a bi-PVA subalgebra of V. According to the general method 
discussed in Section 6.3, L{d) is a pseudodifferential operator of bi-Adler type with 
respect to the bi-PVA A-brackets (7.1) and (7.2). Hence, the local functionals 
fho = 0, fh n = f Res z L n (z), n > 1, are in involution with respect to both 
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PVA A-brackets {• a - }o and {• a -}i, and they satisfy the Lenard-Magri recurrence 
relation 

{Jh n ,u} 0 = {Jh n+1 ,u}i , n>0 ,ueVi. (7.4) 

We thus get the corresponding integrable hierarchy of Hamiltonian equations 
du 

-TT = {fh n ,u}o , n£Z + ,u£V, (7.5) 

CLtn 

which is a bi-Hamiltonian hierarchy when restricted to u £ V±. 

For example, from (7.3) we get the first Hamiltonian functional of the sequence: 


fh 1 = —f Res, L(z) = J -^912921 ■ 


Note that this is a Casimir functional for the first Poisson structure on Vi but not 
on the whole V (cf. Remark 6.9). By a straightforward computation we get 

L 2 (d) = (d 2 + 2 qud + q' n + q\ x - 2 q 12 q 21 - {q 21 +521(911-922)) (<9 + 922) _1 9i2 

— 921(9 + 922) 1 ( — 9i2 + 912(911 ~ 922)) + (92i(9 + 922) 1 9i2) 
and taking its residue we get the second Hamiltonian functional of the sequence: 

1 f 1 

Jh 2 = — i ( z ) = / ^2 (“912921 + 912921(911 _ 922)) • 

The corresponding first two equations of the hierarchy are 

“77~ = _ (9(2 — 912(911 — 922)) , 

at 1 s 

= _ (921 + 921(911 — 922)) , 

at 1 s 

-Jjp 1 = - (912 — 29(2(911 — 922) — 912(911 — 922) + 912(911 — 922) 2 — 2912921) , 

- jp - = ~ (—921 “ 29^(911 — 922) — 921(9(1 — 9(2) “ 921(911 — 922 ) 2 + 2912921) > 
ato s 


^ = ^f= 0 ,forn = 1 , 2 . 

Note that the Lenard-Magri recurrence condition (7.4) holds for u £ V 1 but it 
fails in general for u £ V, cf. Remark 6 . 6 . For example, for n = 0 the LHS of 
(7.4) is zero for every u, while the RHS of (7.4) is not: for u = q 2 \ it is equal to 
921 . On the other hand, one can check that that this hierarchy coincides with the 
homogeneous Drinfeld-Sokolov hierarchy for g ( 2 associated to S = diag(s,0) and 
the diagonal matrix a = diag(l, 0), cf. e.g. [DSKV13]. In particular, this failure of 
the Lenard Magri scheme on V only occurs for n = 0. 


7.2. Second case: S nilpotent. Consider the case when S ' 4 


0 -1 

0 0 


. As 


before, the family of Adler type operators A(d) + eS l , e £ F\{0}, defines the 
compatible PVA A-brackets {-a-}oi given by (7.1), and {• a-} i on V, whose only 
non-zero A-brackets among generators are 


{9ha92i}i — {9i2a922j-i — 1 , {921a9h}i — {922a92i}i — —1 ■ (7.6) 
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Let us consider the canonical factorization S t = IJ where I = and J = 

(0 — l), and consider the corresponding (JJ)-quasideterminant of A(d)\ 


L{d) = \A(d)\ u = (d + q 11 )q 12 1 (d + q 22 ) - 921 


— d 2 + \ 

( 9 ll + 922 


/ 922 \ 

' | 911922 

912 

t <712 

1 2 12 J 

' 9 l 2 ' 

912 


This is a scalar differential operator with coefficients in a differential algebra V 
containing the inverse of qi 2 - 

Since L{d ) is a differential operator, the residue of any its power is zero. On the 
oher hand, assuming that the differential algebra V contains the square root of the 
element q\ 2 , we can consider the square root B{d ) £ V((9 -1 )) of L{d ), and use the 
method described in Section 6.3 to construct a bi-Hamiltonian integrable hierarchy 
for the two Poisson structures associated to the PVA A-brackets (7.1) and (7.6). 

The first few coefficients of the operator B(d) = bid + 6 q + 6-i<9 _1 + ■ ■ • are 


= 


1 


b -1 = - 


v/912 ’ 
1 


bn = 


1 


2Vl 


921 -/ 9 l 2 ~ 


2 V / 9l2 
1 


(911 + 922) — 


9 l 2 


81/912 


+ 


1 


8 V ^ 


(2(911 - 922)912 + 912) - 


4\Z9i2 1 
((911 — Q22) 2 + 2 (g' n — <722)) 
7 


32\/g®2 


( 9 / 


In particular, the first conserved density is hi s = — 2Res z B(z) = — 2b-i and the 
first non-trivial equation of the hierarchy, associated to the corresponding Hamil¬ 
tonian functional, is 


dqn _ d(qn + q 2 2 ) 
dti t B dtiB 

d(qn ~ 922) 


= 0 


,, - 9211/912 + . r— ((911 - 922) 2 + 2 (q'n - ^22)) 

uci.B 4^/512 

(912 + 2(911 - 922)912) + 7 “ — F=?=( 9 i 2) 2 1 


4 V / 9?^ 

d92i 1 


16^ 


dti^B 2 ^ 
1 


(921 — 921(911 — 922)) 


8\/9i2 

1 


( 29 i 292 i — (911 _ Q22) 3 + 2(<7ii — 922) — 912) 


+ 


( 3 ( 9 u ^ 92 2 ) 2 9( 2 - 6(9(1 - 922)9(2 - 2(911 - 922)9(2) 

‘A 

((9ll - 922)( 9 ( 2 ) 2 + 4 9(29i2) - — nr (g'i 2 ) 3 ■ 


32 1 / 9^2 


64 \/<& 


Note that, if we Dirac reduce (cf. [DSKV14]) by the constraint <712 = 1 and qn + 
<222 = 0 and set w = —(<721 + 5(911 — 922 )' + 3(911 — 922 ) 2 ), the above equation gives 
d t™ B = w'. We omit the, rather cumbersome, computation of the next equation 
of the hierarchy, but it is possible to check that, under the same Dirac reduction 
and substitution, we get the KdV equation = \w"' + 2ww'. This is not a 

coincidence: as we will show in [DSKV15new], the Drinfeld-Sokolov hierarchy for 
the affine W-algebra W(gtassociated to a nilpotent element / £ gl N , can 
be obtained, via a limiting procedure, by Dirac reduction of the bi-Hamiltonian 
hierarchy for the affine bi-PVA V/(g( N ). 
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8. Dispersionless Adler type Laurent series and the corresponding 

DISPERSIONLESS INTEGRABLE HIERARCHIES 


8.1. Quasi-classical limit of the algebra of pseudodifferential operators. 

Consider the algebra V((9 -1 )) of pseudodifferential operators over the differential 
algebra V. Recall that the product o on V((d -1 )) is given, in terms of the cor¬ 
responding symbols, by (2.1). Let us consider the family of associative algebras 
14((z -1 )), H G F, whose underlying vector space is V((z -1 )), and the associative 
product is defined by the following deformation of ( 2 . 1 ): 

{A o h B)(z) = A(z + Hd)B(z ), A(z), B(z) G V(( 2 -1 )). 


We denote by [•, -}n the commutator of this associative product. Recall that (cf. 
[LM79]), the quasi-classical limit V qc ((z -1 )) of the family of associative algebras 
{Vft((z _1 ))}?ieF is the Poisson algebra defined as the commutative associative alge¬ 
bra V((z -1 )), with the usual product of Laurent series, endowed with the Poisson 
bracket 


dh 


{A(z),B(z)} qc =—[A(z),B(z)\ h = (d z A(z))dB(z) - (d z B{z))dA{z ). ( 8 . 1 ) 


h- o 


8.2. Dispersionless limit of a family of A-brackets and of an Adler type 
operator. Suppose that A(d) G MatMx n V((9 -1 )) is a matrix pseudodifferential 
operator of H-Adler type , meaning that, for each H G F there exists a A-bracket 
{• a -}r on V such that 


{Aij{z)\A hk {w)}n = A h j(w + HX + Hd)i z (z-w-hX-hd) x | x=g A ik (z 
- A hj (z)i z (z -w-HX- Hd)~ 1 A ik (w). 


hx — hx) 


( 8 . 2 ) 


We will be mainly interested in the scalar case M = N = 1, when the h -Adler 
identity (8.2) for the Laurent series A(z) G V(( 2 -1 )) reads 

{A(z)\A{w)}n = A(w + hX + hd)i z {z — w — HX — hd)~ 1 \ x=d A{z — HX — Hx) 

— A(z)i z (z — w — HX — hd)~ l A{w). 


(8.3) 


Note that (8.2) with H = 1 says that A(d) is of Adler type with respect to the 
A-bracket {• > -}h=i- Note also that the RHS of (8.3) is 0 for h = 0. We define the 
dispersionless limit of the family of A-brackets {• a - }?i as 

{• A‘}disp = A - }?i| K=0 ■ (8-4) 

Then, taking the derivative of both sides of (8.2) with respect to h and letting 
h = 0, we get that A(z) satisfies the following dispersionless Adler identity with 
respect to the A-bracket {• a -}disp: 

{Aij(z)\Ah k )w)) ( \\xp = l z {z w) {d w d z )i^A^j (w) (A -1- d) Ai k (z)) 

+ L z (z - w)~ 2 (A h j(w)(X + d)A ik (z) - A h j{z){X + d)A ik (w)) . 

In the scalar case M = N = 1, the dispersionless Adler identity for a Laurent series 
A(z) G V(( 2 -1 )) reads 

{A(z)AA(w)}dis P = i z {z - w) _1 (5 w - d z )(A(w)(X + d)A(z )) 

+ t z {z - w)~ 2 (A(w)dA(z) - A(z)dA(w)) . 


( 8 . 6 ) 
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8.3. Dispersionless Adler type Laurent series and PVA’s. 

Definition 8.1. A Laurent series A(z) £ Mat^x^ V((z -1 )) is called of dispersion¬ 
less Adler type with respect to a A-bracket {• a -}disp on the differential algebra V 
if equation (8.5) holds. Moreover, we say that A(z) is of dispersionless bi-Adler 
type with respect to a pencil of A-brackets {• a ’}e = {• a - }o + e{- a "}i on V, if 
A(d) + eljv is of dispersionless Adler type with respect to {• a ’}e for every e £ F, 
or, equivalently, if A(z) is of dispersionless Adler type with respect to {• a ’}oj an d 
one has 

{Aij(z)\A hk (w)}i = l z (z - w)~ 1 (6 ik d w A hj (w)X - S hj (X + d)d z A ik (z)) 

+ l z (z - w)~ 2 (5 ik {A hj {w) - A h j(z))\ - 5 h j{\ + d)(A ik (w) - A ik {z))) . 

Remark 8.2. By the arguments in Section 8.2, if A(d) £ MatMxiv V((9 -1 )) is an 
operator of H -Adler type with respect to a family of A-brackets on V, in the sense 
of (8.2), then its symbol A(z) £ MatjvfxJV V((z -1 )) is of dispersionless Adler type 
with respect to the dispersionless limit A-bracket (8.4). In fact, for a scalar Laurent 
series A{z) £ V(( 2 -1 )) the dispersionless Adler identity (8.6) with respect to a A- 
bracket {• a -}disp on V is equivalent to the condition that the S-Adler identity (8.3) 
holds modulo h 2 , with respect to the A-bracket {• a ■}n = h{- a -}disp on V, extended 
by F[[fi]]-linearity to the algebra of formal power series V[[7i]]. (Here we are using 
the fact that, for a scalar operator, the RHS of (8.3) is zero for H = 0.) 

The relation between dispersionless Adler type (scalar) Laurent series and Pois¬ 
son vertex algebras is given by the following result, which is analogous to Theorem 
3.6. 

Theorem 8.3. Let A(z) £ V((z~ 3 )) be a Laurent series of dispersionless Adler 
type, with respect to a X-bracket {• a •} on V, and let Vi d V be the differential 
subalgebra generated by the coefficients of A(z). Then {■ \ •} is a PVA X-bracket on 

Vi. 

If, moreover, A(z) £ V((z -1 )) is of dispersionless bi-Adler type, with respect to 
a pencil of X-brackets {• a •} e = {• a - }o + £ {’ a -} i , e £ F, on V, then {• a -} e , e £ F 
is a pencil of PVA X-brackets on Vi- 

Proof. Of course one can prove Theorem 8.3 by a straightforward but long compu¬ 
tation. Instead, we shall use a different argument. The same proof of Theorem 3.6 
(cf. the proofs of [DSKV15a, Lem.2.2 & Lem.2.5]) can be adapted to prove that, 
if A{9) £ V((9 -1 )) satisfies the fi-Adler identity (8.3) with respect to a family of 
A-brackets {■ \-}n, ft £¥, then the skewsymmetry and Jacobi identity axioms hold 
for these A-brackets for every H £ F. Moreover, still the same proof shows that, if 
the A-bracket {• a -}fi. has values in fiV[[ft]] and A(d) £ V((9 -1 )) satisfies the ?i-Adler 
identity (8.3) modulo H 2 , then the skewsymmetry axiom holds modulo h 2 : 

{A(z)xA(w)} h + {A(w)-x-dA(z)} h = 0 mod ( h 2 ), 
and the Jacobi identity holds modulo h 3 : 

{A(zi)\{A(z2) k ,A(z3)}fi}h ~ {A(z2)i_i.{A(zi)xA(z3)}fi}h 
- {{A(zi)AA(z 2 )}fi A+M A( 23 )} fi = 0 mod (e 3 ). 

The claim follows from Remark 8.2. The arguments for the second statement are 
the same. □ 

Example 8.4. Consider the symbol of the generic matrix pseudodifferential (resp. 
differential) operators from Example 3.5: 

M M 

Lmn{z)= ^2 UjZ 3 (resp. L( MN) (z) = ^17^) with U M = Ijv , (8.8) 

j =—oo j=0 
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/ \N 

where Uj = {Uj- a b) € Mat^xArV for every j < M, and V is the algebra of 
differential polynomials in the infinitely many variables zij ;a b, where 1 < a, b < N 
and — oo < j < M (resp. in the finitely many variables Uj ;a b , where 1 < a, b < N 
and 0 < j < M). For the same reasons that Lmn{9) (resp. L^ MN ^(d)) is of bi-Adler 
type, it is immediate that, for every e € F, Lmn(z) + etw (resp. L( MN )(z) + eljv) 
is a Laurent series (resp. polynomial) of h -Adler type with respect to a pencil of A- 
brackets {• a •} h,e = {• a ‘jfi.o + e{ - a e € F, depending on the parameter H £ F, 
on V. These A-brackets are (uniquely and well) defined by the 7i-Adler identity: 

{LMN-,ij{z)\LMN-hk(,w)}h,e 

= ( LMN-hj(,w+hX+Hd) + eShj)bz{z-w-hX-Hdy 1 \ x=g (L M N-,ik{z-hX-hx) + eSik) 
- ( L MN -hj{z) + e6 h j)i z (z -w-h\- hd)-\L M N-,ikW + e< 5*fc) > 

(8.9) 


(and the same equations for L( mn) in place of Lmn)- It follows by the arguments 
in Section 8.2 that, taking the dispersionless limit (8.4), the Laurent series Lmn(z) 
(resp. the polynomial L^ M n){z)) is of dispersionless bi-Adler type with respect to 
the pencil of A-brackets {• a -}disp,e = - a '}n,e\ h=0 - By Theorem 8.3 we have that, 

in the scalar case N = 1, this is a pencil of PVA A-brackets on V. The two PVA 
structures are given, respectively, by (denoting L(z) either Lm\{z) or 

{L(2:)AL(ui)}disp,o = ^{z - w)~ 2 (L(w)dL(z) - L{z)dL{w)) 

+ i z (z - w)~ 1 (d w - d z )(L(w )(A + d)L(z)) , 

and 

{L(2)AL(w)}di sp ,i = l z {z - w)~ 2 d(L(z) - L(w)) 

+ i z (z - w) -1 (d w L(u j)A - (A + d)d z L(z)) . 

We can rewrite the two A-brackets (8.10) and (8.11) in terms of the generators 
In the Laurent series case we have L(z) = Lmi(z) = Y2j=-oo u j z '’ 
(um = 1), and the two brackets (8.10) and (8.11) are (—oo < i, j < M — 1) 

{uiA%}disp,o = (j + l)tij+i(A + d)u i+ i 

M—i—2 M—i—2 


( 8 . 10 ) 


( 8 . 11 ) 


+ U-*) J2 Uj~e(\ + d)iii+g +2 ~ (2A -I- d ) ^ ' (l + l)v,j—£Ui+£+2 ■ 


( 8 . 12 ) 


e=o 


1=0 


and 

( -((f+l)(A-t-9)-|-(j-|-l)A)wj + j + 2 for i,j > 0 s.t. i+j+2 < M 
{u ix Uj}di S p,i = < ((«+l)(A+<9) + (j + l)A)u i _ h7 +2 for i,j < -1 
[ 0 otherwise. 

(8.13) 

The Poisson structure (8.13) for M = 1 was found in [KM78]. In the polynomial 
case we have L(z) = L(mi){z) = Yhj=o u j z j '> and the t wo brackets (8.10) and (8.11) 
are (0 <i,j < M — 1) 

{MiA M i}disp,0 — (j 4“ l)«j+l(A -{- d^Uid-l 

min{M—i—2,j} min{M—i—2,j} 

+ u - *) 5Z Uj—t {A + d)ui+i+2 ~ (2A + d) 'y ' (£ + l)iij-eUi+e+2 , 

t =o i =o 

(8.14) 

and 

{«.«L„P,1 = { -((ViKA+ai + O+LA),,,^ (8.,5) 
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Note that this bi-PVA is a quotient of the previous one by the differential ideal 
generated by {uj | j < 0}. 

In the matrix case TV > 2 it is not true, in general, that {• a -}disp,e is a pencil of 
PVA structures on V, since the limit H —> 0 of the RHS of (8.9) is, in general, not 
zero. It is true, however, in the polynomial case of order 1, i.e. for L^ 1N )(z). In this 
case, the 0-th A-bracket {• a -}disp,o is a PVA A-bracket, which is given explicitly, 
letting L (1N) {z) = 6ijZ + Uij, by 

\Uhk}disp,0 = dhjdikX ; (8.16) 

while the 1-st A-bracket {• a ‘}disp,i is identically zero. 

8.4. Integrable hierarchies associated to dispersionless Adler type Lau¬ 
rent series. In this section we show that, given a Laurent series of dispersionless 
Adler type with respect to a PVA A-bracket, there is an associated sequence of 
Hamiltonian functionals in involution, namely the dipersionless analogous to Theo¬ 
rem 5.1, and that given a Laurent series of dispersionless bi-Adler type with respect 
to a pencil of PVA A-brackets, the corresponding sequence of Hamiltonian function¬ 
als satisfies the Lenard-Magri recurrence relation with respect to the given Poisson 
structures, namely the dipersionless analogous of Theorem 6.4. 


Theorem 8.5. Let V be a differential algebra with a PVA X-bracket {-a-}- Let 
A[z) £ V((t; _1 )) be a Laurent series of degree K > 1 with invertible leading coeffi¬ 
cient, and let B(z) £ V((;r -1 )) be its K-th root, i.e. A(z) = B(z) K . Assume that 
A(z) is of dispersionless Adler type with respect to the X-bracket {• a •}• Define the 
elements h n £ V, n £ Z, by 

h n = —Res z B(z) n for n / 0 , h 0 = 0. (8-17) 

M 


Then: 


(a) All the elements Jh n are Hamiltonian functionals in involution: 

{ Jh m , fhnj = 0 for all m,n £ Z. 


(8.18) 


(b) The corresponding compatible hierarchy of Hamiltonian equations is 

^ = {fh n , A(z)} = {(. B(z) n )+, A{z)} qc ,n£ Z (8.19) 

(the bracket {•, -} qc in the RHS is given by (8.1),), and the Hamiltonian func¬ 
tionals Jh n , n £ Z + , are integrals of motion of all these equations. 

(c) If, moreover, A(z) is a Laurent series of dispersionless bi-Adler type with respect 
to two PVA X-brackets {• a - }o and {• a - }i on V, then, the elements h n , n £ Z+, 
satisfy the Lenard-Magri recurrence relation: 


{fh n ,A(z)} 0 ={fh n+K ,A(z)} 1 ={(B(z) n ) + ,A(z)} qc , n£ Z. (8.20) 
Hence, (8.19) is a compatible hierarchy of bi-Hamiltonian equations. 


The proof of Theorem 8.5 relies on the next Lemmas 8.6 and 8.7. Lemma 8.6 
is analogous to Lemma 5.5 and it is an obvious consequence of the Leibniz rules 
(recall that the algebra V((z -1 )) is commutative). Lemma 8.7 is straightforward. 
We state them for completeness. 

Lemma 8.6. For a £ V, n £ Z, and A(z) and B(z ) £ V(( 2 -1 )) as in Theorem 
8.5, we have 

{*nA«}|A =0 = -R- e Sz{^(2)xa}(| x=a B(^)" _if ) , 

/{aAM | A=0 = ~ J Res z B(z) n ~ K {a x A(z)}\ x ^ 0 


( 8 . 21 ) 
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Lemma 8.7. For A(z),B(z) £ V{{z 1 )), we have 

f Res,{A(.z), B(z)} qc = 0. (8.22) 

Proof of Theorem 8.5. Applying the second equation in (8.21) first, and then the 
first equation in (8.21), we get 

{fh m ,fh n } = J Res z Res w B(w) n ~ K {A(z) x A(w)}{\ x=g B(z) m ~ K ). (8.23) 

We can now use the dispersionless Adler identity (8.6) and the fact that A(z) = 
B(z) k to rewrite the RHS of (8.23) as 


J Res, Res™ t z (z — w) 2 (B(w) n B(z) m K dA(z) — B(z) m B(w) n K dA(w )) 

+ J Res z Res w (d(B(z) m )B(w) n - K d w A{w) - B(w) n d(B(z) m - K d z A(z))) . 

(8.24) 


Note that 


K 


B(z) m ~ dA(z) = —dB(z) m , 
m 


(8.25) 


and similarly with d replaced by d z or d w . Hence, we can rewrite (8.24) as 

/ Res z Res™ l z (z — w)~ 2 (—B(w) n dB(z) m — —B(z) m dB(w) n \ 

J \m n J 

+ f Res 2 Res™ l z (z — w)~ x (-d{B(z) m )d w B(w) n --B(w) n d{d z B(z) m )] . 
J \n m J 


(8.26) 


By the fact that d and d z commute and using integration by parts, we get 

Res z l z (z — w)~ l d(d z B{z) m ) = Res 2 i z (z — w)~ 2 dB{z) m . (8.27) 

Furthermore, by (2.6), we have 

Res, l z (z - w)~ 1 dB(z) m = {dB{w) m ) + = d(B(w) m )+ , (8.28) 

where in the last identity we used the fact that d acts componentwise on elements 
of V((w -1 )), and 

Res, t z (z - w)~ 2 B(z) m = - Res, (d~t z (z - w)” 1 ) B(z) m 

= (d w B(w) m ) + = d w (B(w) m )+ . [ "' 2U) 

Here, in the second identity we used integration by parts and equation (2.6), and 
in the last identity the fact that taking the positive part of elements in V((u; _1 )) 
commutes with the partial derivative d w . Hence, using equations (8.27), (8.28) and 
(8.29) we can rewrite equation (8.26) as 

- [ Res w {d(B{w) m ) + d w B(w) n - d w (B{w) m ) + dB{w) n ) 

n J 

which is zero by Lemma 8.7 thus proving part (a). 

For part (b), by the first equation in Lemma 8.21 we have 

{Jh n ,A(w)} = {h nA A(w)}| A=0 = - Res z {A(z) x A(w)}(\ x=d B(z) n ~ K ) . (8.30) 
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By the dispersionless Adler assumption (8.6), the identity A = B K 
equation in (8.25), the RHS of (8.30) becomes 

- Res z l z (z - w)~ 2 A(w)dB{z) n - B{z) n dA(w)^ 

- Res z t z (z - w) _1 ^ d w (A(w))dB(z) n - A(w)d (d z B(z) n )^ 
= d w (B(w) n ) + dA(w) - d w (A(w))d(B(w) n ) + . 


and the first 


(8.31) 


In the second equality we used equations (8.27), (8.28) and (8.29). This proves 
(8.19) and completes the proof of part (b). 

Finally, we prove part (c). By the first equation in (8.21) we have 

{fh n ,A(w )}i = {h n \A{w)}i\ x=0 = -Res 2 {A( 2 ;) x A('U))}i(| a , =a -B” _ ' R ’(^)) ■ (8.32) 

Applying the dispersionless bi-Adler identity (8.7) and (8.25) we can rewrite the 
RHS of (8.32) as 


K 


- Res z t z (z - w)~ l — dB(z) n - B(z) n ~*dA(w) 
n 


— Res z ^ d w A(w)dB(z) n K 
= {(B(z) n - K )+,A(w)} qc . 


K 


dd z B(z) r 


(8.33) 


In the first equality we used again equations (8.27), (8.28) and (8.29). Comparing 
(8.19) and (8.33), we get (8.20). □ 


Remark 8.8. Note that the Poisson structures associated to the generic Adler type 
series in Example 8.4 coincide with the corresponding dispersionless limit in [DZ01] . 
Hence, the hierarchy (8.19) can be obtained as the dispersionless limit, as defined 
by Dubrovin and Zhang, of the hierarchy (5.5). 

8.5. Examples. 

Example 8.9 (dKP hierarchy). Let V be the algebra of differential polynomials in 
the variables {uj | i < 0}. The Laurent series L(z ) = z + Y^i<o Uizi °f dispersion¬ 
less bi-Adler type with bi-PVA structure defined on generators by equations (8.12) 
and (8.13) for M = 1, see Example 8.4. We can compute the first few integrals of 
motion Jhk, k> 1, directly from the definition (8.17): 

fh 1 =-fit- 1 , fh 2 = -J(u -2 + u 0 u-i), 

f ll 3=-f (u_3 + 2u 0 U-2 + U 2 _ 1 + UqU_ 1 ) , . . . 

To find the corresponding bi-Hamiltonian equations, we use Theorem 8.5. We have 
L(z)+ = z + ito, L 2 (z)+ = z 2 + 2 uqz + 2m_i + Uq, L 3 (z)+ = z 3 + 3uoz 2 + 3(u_i + 
iiy)z + ( 3u- 2 + 6iio‘U-i + Uq). Hence, 

= dL(z) - d z L(z)u' 0 , 

- = (2~ + u 0 )dL{z) - 2d z L(z){u' 0 z + u'^ + u 0 u ' 0 ), ^ 

= (3 z 2 + 6uoz + 3(u_i + Uq)) dL(z) 
ats 

- d z L(z ) (3 u' Q z 2 + 3 (u'_i + 2u 0 u' 0 )z + 9(3 u_ 2 + 6u 0 u_i + Uq)) ,... 
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Remark 8.10. Since u o is central for the A-bracket (8.13), we have ijj 2 - = 0, for 
n £ Z + . Hence, we may assume u o = 0 in (8.34). The second equation in (8.34) 
then becomes 

= 2(«*,_]. - (fe + l)uk+iu '_ x ), fc<-l, (8.35) 

which is the Benney equation for moments, see [KM78, LM79|. Moreover, consider 
the first two equations in the second system of the hierarchy (8.34), and the first 
equation in the third system of (8.34). After eliminating the variables U -2 and u -3 
and relabeling t -2 = y, £3 = t and u = 2 zt_i, we get 

3 u V y = (4 ut — 6uu ')', (8.36) 

which is known as the dispersionless Kadomtsev-Petviashvili equation. 

Remark 8.11. It has been shown by Radul [Rad87] that there are infinitely many 
bi-Poisson structures for the KP hierarchy, see also [DSKV15a, Rem.3.14]. In fact, 
we have infinitely many bi-Poisson structures for the dKP hierarchy as well, corre¬ 
sponding to the bi-PVA structures defined by (8.12) and (8.13) on the algebra of 
differential polynomials in the variables {zq | i < M — 1}, for every M > 1. These 
bi-PVA structures are obtained as the dispersionless limit of the bi-PVA structures 
for the KP equations. 

Example 8.12 (The dispersionless M-th KdV hierarchy). Let us consider the poly¬ 
nomial L(mi)(~) = 0 u i z% 1 w bh mm = 1. It is of dispersionless bi-Adler type for 

the A-brackets (8.14) and (8.15) defined on the algebra of differential polynomials in 
the variables {zq | 0 < i < M— 1}, see Example 8.4. It is proved in [DSKV15a] that, 
after setting um- 1 = 0 (which is possible since um-i is central for the A-bracket 
(8.15)), the integrable bi-Hamiltonian hierarchy ( 6 . 6 ), associated to the Adler type 
differential operator Lt M is the M-th KdV hierarchy. Hence, the integrable 
bi-Hamiltonian hierarchy (8.20) is the dispersionless M-th KdV hierarchy. 

9. Adler type pseudodifferential operators for double Poisson 

VERTEX ALGEBRAS 

In this section, by a differential algebra V we mean a (possibly noncommutative) 
unital associative algebra with a derivation d. We refer to [DSKV15b] for the 
definition and the main properties of a structure of a double Poisson vertex algebra 
on V. 

Given a differential algebra V, we define for each positive integer N a unital 
commutative associative differential algebra Vat, generated by elements a^, a £ V, 
i,j £ {1,..., N}, subject to the following relations (k £ F, a, b £ V, i, j = 1,..., N ) 

N 

( ka'jjj — ka.ij , ( 0 . T tyij a^j - 1 - b^j , (zz^)^ — ^ ) O'ikbkj > (^*-0 

fc—1 

with derivation d : Vat —> Vn given by daij = ( da)ij. Recall from [DSKV15b, 
Prop.3.20, Thm.3.22] that, if V is a double Poisson vertex algebra, with 2-fold 
A-bracket ({• a •§, written (in Sweedler’s notation) as 

ia x bj = J2 ( a ^)' ® (onb)" X n , (9.2) 

n(EZ-|_ 

then we have a PVA structure on Vat, N > 1, with the A-bracket given by 

{a.ij\bhk}= ( a n b )hj( a nb)ik A n . 
nGZ-|_ 


(9.3) 
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Recall also that, given a double PVA V, we have a well defined Lie algebra 
structure on the quotient space V/([V, V] + dV) defined by 


U a J b } = f m 


a=o = f(aobY(a 0 by 


a,b G V, 


(9.4) 


where f : V —> V/([V, V] + dV) denotes the projection map and m : V (g) V —> V 
denotes the multiplication map. Furthermore, we have a representation of the Lie 
algebra V/([V, V] + 9V) on V, with the Lie algebra action given by 


{/a, b} = m 


= (a 0 &) , (a 0 6)'' 


A=0 


a, b G V . 


(9.5) 


We shall also denote (for clarity of the exposition), for every N > 1, J N : Vn —> 


V/v/9V/v the canonical quotient map and by {• 
(9.3) on the quotient space Vn/RVn- 


the Lie bracket induced by 


Proposition 9.1. (a) There is a well defined Lie algebra homomorphism ipN 


V/([V, V] +9V) -> V N /dV N given by VN(fa) = f N J2i=i a u- 
(b) For every i,j = 1,..., N, the map a <—> aij, from the module V over the Lie 
algebra V/([V, V] +9V) to the moduleVN over the Lie algebra Vn /9Vn , is com¬ 
patible with the Lie algebra homomorphism pn ■ V/([V, V] + dV) —> Vn/9Vn 
defined in (a), meaning that 


{fa,b}ij = {ip N (Ja), bij} 
Proof. Straightforward. 


for all a,b £ V . 


(9.6) 

□ 


Remark 9.2. In fact, if V is a double Poisson vertex algebra, then we have a Lie 
conformal algebra A-bracket on V/[V, V], given by {tr(a)A tr(5)} = tr (m^aA^])), 
a, b G V, where tr : V — > V/[V, V] is the canonical quotient map, and we have a Lie 
conformal algebra homomorphism V/[V, V] —> Vn, given by tr(a) i—> y )^!_ 1 an. 

Definition 9.3 ([DSKV15b]). Let V be a differential algebra endowed with a 2-fold 
A-bracket a'S- We call a pseudodifferential operator A(9) G V((9 -1 )) of Adler 
type for ({• a •]) if the following identity holds in V® 2 [A,^]((^ _ \w !)): 


§A(2)aA('ui)3- = A(w + A + d) ® i z (z — w — X — d) 1 A*(—z + A) 
— A(z) (g) i z (z — w — A — 9) _1 A(u;). 


(9.7) 


Furthermore, we say that A(d) is of bi-Adler type with respect to the two 2-fold 
A-brackets -f-A'S'O and j[‘A‘]}i if it is of Adler type with respect to the 2-fold 
A-bracket ({• a -]j-o> Le. (9.7) holds, and 

l z (z — w — A — 9) _1 (A*(— z + A) — A(w)) 

_i /a / a/ ^\ (9.8) 


^A(z) a A(u;)Ji = 1< 

+ i z (z — w — A) - 


(A(w + A) — A(z)) ® 1. 


The following results are a generalization of the results in [DSKV15b, Sec.5] and 
can be proved adapting the proofs of Theorems 3.6, 5.1, 6.3 and 6.4 to the double 
Poisson vertex algebra setting. 

Theorem 9.4. Let V be a differential algebra, endowed with a 2-fold X-bracket 
j{-A‘Jj'. Let A(d) £ V((9 -1 )) be a pseudodifferential operator of Adler type with 
respect to ({• a •})■ Let Vi C V be the differential subalgebra generated by the coeffi¬ 
cients of A(d). Then: 

(a) V± is a double Poisson vertex algebra (with double X-bracket ({• a •]},)■ 

(b) For B(d) G V((9 -1 )) a K-th root of A (i.e. A(d) = B{d) K for K G Z\{0}j 
define the elements h n ^s G V, n G Z, by 

h n ,B = q—r Res z B n (z ) for n/0, h 0 = 0 . 


(9.9) 
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Then all the elements fh n ^B are Hamiltonian functionals in involution: 

{/ f h nt c} = 0 , for all m,n £ Z, B,C roots of A . 

The corresponding compatible hierarchy of Hamiltonian equations is 

-T7 — = {f h n,B,A(z)} = {(B n ) + ,A](z), n £ Z, B root of A, (9.10) 

atn,B 

and the Hamiltonian functionals jh Ut c, n £ Z+, C root of A, are integrals of 
motion of all these equations. 

Furthermore, let us assume that A(d) £ V((<9 -1 )) is of bi-Adler type with respect to 
the 2-fold X-brackets ({• a - Jo and S' a •}} i on V. Then: 

(c) V is a bi-double PVA with the 2-fold X-brackets ({• a - Jo an d a ■ Ji- 

(d) For a K-t.h root B(d) £ V((<9 -1 )) of A(d), the elements h U} B £ V, n £ Z + , 
given by (9.9) satisfy the generalized Lenard-Magri recurrence relation (6.6). 
Hence, (9.10) is a compatible hierarchy of bi-Hamiltonian equations. 

Note that by Proposition 9.1(a), the sequence h U) B £ V, n £ Z + , defined by 
(9.9) gives rise to a sequence 

h n,B = Tn{J h n , B ) eVjv, n £ Z+, N > 1, 

of Hamiltonian functionals in involution, and by Proposition 9.1(b) we get the inte¬ 
grate Hamiltonian hierarchy in Vn corresponding to (9.10), which can be regarded 
as the Hamiltonian hierarchy for the entries of the matrix hierarchy (9.10). 

Example 9.5 ([DSKV15b, Ex.3.23]). Consider the noncommutative algebra of 
differential polynomials in two (noncommutative) variables u and v endowed with 
the following 2-fold PVA A-bracket: 

]{uam]} = 1 ® u — ii <8> 1 + c(l (g> 1)A , ceF, v central . 

The elements ho = 1, h n = ^(u + v) n for n > 0, are in involution, and the corre¬ 
sponding integrate hierarchy of Hamiltonian equations is: 

yp - = v(u + v) n U — u{u + v) n v + cd(u + 1 ))" +1 , yy^- = 0 , U £ Z + . 

CLin din 

For every N > 1, Vat is the commutative algebra of differential polynomials in the 
variables Uij and Vij, 1 < i, j < N, endowed with the following PVA A-bracket 

{uijyUhk} = ShjUik - 5 ik Uhj + cSikShjX , central for every i,j . 

The corresponding Hamiltonian densities are 

T N 71 / \ 

h$ = N, C = J(Ai(Ao n> 1, 

i,j =1 fc=0 ' 


where in the RHS above we have used equation (9.1). The Hamiltonian equations 
become 


dui 

dt n 

dui 

dt r 


N 


= cd({u + v) n+1 )ij + ^2 ( Vih((u + v) n )hk'Ukj - U ih ((u + v) n ) hk v k j) , 
h,k= 1 


= 0, 1 < i, j < N, n £ Z 


+ ■ 
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Appendix A. Square roots of <9l 2 + Q G Mat 2x2 V((<9 -1 )) 

Let K be a differential field with subfield of constant F. Let Q = (fiji) i ._ 1 G 
Mat 2x2 /C. The matrix differential operator A(d) = 91 2 -f Q G Mat 2x2 /C[<9] has 
Dieudonne determinant of degree 2, and it is natural to ask whether A{d) has a 
square root in Mat 2x2 /C((9 -1 )) or over an extension of 1C. 

Lemma A.l. In a differential field extension K. of 1C, there exists an invertible 
matrix B G MatTVxAf £ such that 

B' + QB = 0. (A.l) 

Proof. Equation (A.l) can be rewritten as the following systems of TV 2 1-st order 
linear differential equations in the TV 2 entries bij of the matrix B: 

N 

b 'ij + = 0 > 

1=1 

which has an TV 2 dimensional space of solutions in a differential field extension K. 
of 1C. □ 


Lemma A. 2. All the possible square roots of the matrix l 2 d are the matrix pseu¬ 
dodifferential operators with constant coefficients A G Mat 2X 2 F((9 -1 )) such that 
tr(A) = 0 and det(A) = d. Namely, they are the matrices of the form 

a(d) /3(d) 

(3(d) 1 (a(d) 2 — d) —a(d) 

with a, (3 G F((9 -1 )) and (3(d) 0. 

a(d) b(d) 


A = 


(A.2) 


Proof. Let A = 


c(d) d(d) 

pseudodifferential operators a, b,c,dC /C((9 -1 )) satisfy 

a(d) 2 + b(d)c(d) = d 
a(d)b(d) + b(d)d(d) = 0 
c(d)a(d) + d(d)c(d) = 0 
c(d)b(d) + d(d) 2 = d. 


, be a square root of d 1 2 . In other words, the 


(A.3) 


It cannot be b(d) = 0, otherwise from the first equation we get a(d) 2 = d, which is 
impossible by obvious degree considerations. Hence, b(d) 0, and from the second 
equation we get 

d(d) = -b(d)~ 1 a(d)b(d). (A.4) 

Substituting (A.4) in the other three equations of (A.3), we get 
a(d) 2 + b(d)c(d) = d 

c(d)a(d) - b(d)~ 1 a(d)b(d)c(d) = 0 (A.5) 

c(d)b(d) + b(d)~ 1 a(d) 2 b(d) = d. 


If we multiply the third equation of (A.5) on the left by b(d) and on the right by 
b(d) _1 , we get 

b(d)c(d) + a(d) 2 = b(d)dHd )- 1 , 
which, when compared to the first equation of (A.5), implies 

b(d)d = db(d ), (A.6) 

namely b(d) G F((9 -1 )). Next, if we multiply on the left by b(d) the second equation 
of (A.5) we get 


b(d)c(d)a(d) = a(d)b(d)c(d ), 
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and this equation can be rewritten, using the first equation of (A. 5), as 

<9a(<9) = a(d)d , 

namely a(<9) £ F((<9 -1 )). To conclude, a(<9) and b(d), being with constant coeffi¬ 
cients, commute, and therefore equation (A.4) gives d( d) = —a(d). Moreover, the 
first equation of (A. 5) gives c(<9) = — b(d)~ 1 (d — a{d) 2 ). This proves (A.2) and the 
claim. □ 

Remark A.3. It is not true, in general, that a matrix pseudodifferential operator 
A(d) £ Mat a r x iv V((9 -1 )) whose degree of the Dieudonne determinant is equal to 

/ d 3 0 \ 

A > 1 admits a A-th root. For example, the matrix I g I has Dieudonne 

determinant of degree 4. But the same computations used in the proof of Lemma 
(A.2) show that this matrix cannot have any square root (equation (A. 6) is replaced, 
in this case, by equation b(d)d 3 = db(d ), which does not admit any solution, by 
obvious degree considerations). 

Proposition A.4. Let B £ Mat 2 X 2 A be any invertible matrix solving (A.l). Then 
all possible square roots of the matrix A(<9) = d + Q are of the form 

Ai(d) = BA{d)B- 1 , (A.7) 

where A(<9) is an arbitrary root of d I 2 , as in Lemma A. 2. 

Proof. Note that, under the assumption (A.l) on A?, we have BdB~ x = d + Q. 
Hence, if A^(d) is a square root of A(d ), then 

(A?{d)) 2 =8 + Q = BdB~ 1 , 

and multiplying on the left by B~ 1 and on the right by B , we get 

(B~ 1 A^(d)B) 2 = dl 2 . 

In other words, B~ 1 A^(d)B is a square root of <91, i.e. it is B~ X A^B = A(9), as 


in Lemma A.2. The claim follows. 


□ 


Remark A.5. Consider the local density functions obtained by taking residue of 
fractional powers of the operator A(d) (cf. (5.3)): 

h n = tr Resg A% (d ), n > 0 , 

where Ai(d) = (<912 + Q)i is one of the square roots given by Proposition A.4. 
For even n = 2k, A% (<9) = Bd k B~ l , so the residue is zero unless n = —2. On the 
other hand, for odd n = 2k + 1, we have 

h n = tr Resa BA n (d)B~ 1 = tr Res^ BA(d)d k B~ 1 . (A.8) 

By the cyclic property J trResaAR = f trResg BA, and since, by Lemma A.2, 
A(<9) has zero trace, the element h n lies in the image of <9 (in 1C). 
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